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Analysis and Standardization of 
Astrodynamic Constants 


M. W. Makemson;’ R. M. L. Baker, Jr.;2? and G. B. Westrom’ 


stract 


\n analysis is presented of heliocentric, geocentric, seleno- 
tric, planetocentric, and atmospheric astrodynamic 
stants. A thorough survey of past data is presented and 
et of adopted values for all constants and their errors is 
posed in order to establish a standard. 


troduction 


Classical astronomical constants, however precisely 
ermined, do not suffice for predicting trajectories 
space vehicles. Motions of close satellites, for ex- 
ple, have brought to light higher order terms in 
» Harth’s gravitational potential which have no 
act on the distant Moon. For successful lunar and 
erplanetary missions more exact values of the solar 
rallax, the astronomical unit in laboratory units, 
1 the distances, diameters, figures, masses, tem- 
atures and atmospheres of the Moon and target 
nets are clearly essential. 

The new astrodynamic constants, associated with 
wee vehicle trajectories, may be viewed in three 
ects: (a) the advance determination of the most 
ybable values; (b) pre-knowledge of the error to be 
yected in the final outcome and how to correct it in 
ht; and (c) improvement of the constants on the 
sis of further observations, usually by a least squares 
ferential correction of the orbit. For example, the 
sent uncertainty in the Earth’s radius may lead 
a 78-mile error on lunar impact as shown by Tross 
0)° and an inaccuracy of 0”01 in the solar parallax 
y produce a divergence of 50,000 miles from Mars. 
An occasional review and comparison of the various 
erminations of astrodynamic constants is clearly a 
sessity, and it is hoped that the adoption of con- 
rent values, as suggested at the end of each table 
the present paper, may serve to standardize compu- 
ions. 


“This work was supported by the United States Air Force, 
part, under Contract AF 49(638)-498 monitored at UCLA 
the AFOSR. Presented at the Seventh Annual Meeting of 
American Astronautical Society January 17, 1961. Manu- 
ipt received September 1960, revised December 1960. 
‘University of California, Los Angeles. 

Air Force Ballistic Missile Division, Directorate of Air 
ce Space Boosters, Los Angeles. 

Aeronutronic Division of Ford Motor Company, Boston. 
Numbers in parentheses refer to references found at the 
elusion of the paper. 


1. Heliocentric Constants 


In a heliocentric orbit the solar parallax, Ig, and 
the gravitational constant, k, , are of first importance, 


II (radians) = = (1) 
® 


where a, is the Earth’s equatorial radius and a@ its 
mean distance from the Sun. The astronomical unit 
was originally defined as the semi major axis of the 
Karth’s orbit and was so used by Gauss, neglecting 
perturbations, in his classic determination of the 
numerical value of k,, the “Gaussian constant.” In 
order to preserve Gauss’ value and eliminate necessary 
revisions of planetary tables as more accurate funda- 
mental constants became available, the definition of 
the astronomical unit was changed to “the mean 
distance of a fictitious unperturbed planet having the 
mass and sidereal period that Gauss adopted for the 
HKarth” when he calculated k, from Kepler’s third 
law, 1.e., 


z PAZ i c 
where m is the mass of a planet in units of the Sun’s 
mass; @ its mean distance from the Sun in astronomical 
units, and P its sidereal period. As Gauss’ value of 
k, is retained for convenience by international agree- 
ment, it is clear that the Earth’s mean distance in 
astronomical units (a.u.) is subject to improvement as 
more exact determinations of the solar parallax and 
the mass and sidereal period of the Earth become avail- 
able. Predicted planet positions, also dependent on k, , 
are thus improvable up to the limits imposed by 
knowledge of the Earth’s mass and period of revolution, 
subject to the errors of observation in an individual 
case. The general limitation is of the order of 8 or 9 
significant figures. The adopted values for these funda- 
mental heliocentric constants are 


k, = 0.017,202,098,95 radians/mean solar day 
= 3,5487187,606,965,1 of arc/mean solar day 

ade@ = 1.000,000,03 a.u., Newcomb (76) 1895; 
American Ephemeris (5). 

Present inaccurate values of the solar parallax 


greatly restrict the accuracy of the ratio of the as- 
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tronomical unit to laboratory units, thus downgrading 
the 8 or 9 figure accuracy of the astronomical system 
to the 4 or 5 figure significance of the ratio. Table 1 
illustrates the situation. In Table 1 Rabe’s value is 
adopted with a larger uncertainty in order to partially 
encompass the values recently determined by Price 
et al. (as recomputed by Westrom and Arsenault 
(40)) and McGuire et al. (will not encompass Spencer 
Jones’ value, however). The error is chosen not only 
on the basis of the probable experimental errors quoted 
by the investigators; but also upon an estimate of 
systematic errors that must underlie the diversity in 
the results. The same philosophy will govern errors 
adopted in subsequent tables. 

Neglecting a small uncertainty in a,, the Harth’s 
equatorial radius, the divergences of interplanetary 
trajectories corresponding to an inaccuracy in the solar 
parallax, expressed as the ratio 


Allo 0.002 
TW — , ee 0.0 93 « : | 
ine ( 8.798 00,23 adopted 19 it) (3) 


are found in Table 2. 


2. Geocentric Constants 

Mathematical terms in the gravitational potential 
representing departures of the reference spheroid (ap- 
proximating the Earth’s figure) from a concentrically 
homogeneous sphere are of two types: zonal har- 


monies, expressing the meridional ellipticity and a 
metry, and tesseral or sectorial harmonics in longi 
which due to the Earth’s rotation must be expr 
as functions of the time. 

The form of the Earth’s gravitational po 
given by Baker and Makemson (7) (Equatio 
page 90) as: 


2 
® eau) Ee on — 3 sin’ 6) 
if oT 


40) (ae aecin aie 
ete 
+! 3 _ gon’ sq 35's 0 eee | 
a0 
4 


should, probably, be modified in the light of Brouw 
suggestion (10) to read ‘ 


kom J» ey . 

sige 14+ —((1—3 sin 3) 
r 2 3 

<= Js (3 — 5 sin’ 6) sin 6 

28 

Ahi 2292, ae 
= += (3 — 30 sin 6 +35 sim’ 6) 

8r4 

Je 


— ©* (15 — 70 sin’ 6 + 63 sin’ 6) sin 6 


* Tere (5 — 105 sin® 6 + 315 sin’ 6 — 231 sin’ 6) + § 


TABLE 1 
The Solar Parallax and Earth’s Mean Distance from Sun 


Solar Parallax® Author Method CO 
(X 106 km 
8.8036 + 07004,6 Gill (85) 1897 triangulation_on Victoria, Iris, Sappho | 149.437 
8.806 + 0.004 Hinks (44) 1904 triangulation on Eros . 149.40 
8.803 + 0.004 Spencer Jones (97) 1928 | radial velocities 149.45 
8.799 + 0.001 Witt (107) 1935 dynamical—Eros (ef. section 3.1) 149.52 
8.790 + 0.001 Spencer Jones (98) 1941 triangulation on Eros 149.67 
8.805 + 0.007 Adams (1) 1941 radial velocities of stars 149.41 
8.798,4 + 0.000,4 E. Rabe (85) 1949 | dynamical—Eros 149.526 
8.792 + 0.003 Brouwer (9) 1950 occultations» 149.63 
8.801,84 + 0.000,05 Price (84) 1959 | radar echoes from Venus (recomputed 149 .467 
by Hilton and Arsenault) 
pee + 0.000,1 MeGuire, Morrison, Wong (69) 1960 Doppler data from Pioneer V* 149.54 
8.797,9 + 0.000,2 de Vaucouleurs (30) 1961 Adopted from mean of several determi- 149.53 
: nations 
8.798 + 0.002 Adopted (1961) 


(Newcomb’s value determined in 1895 of 880 is : : derioan acim pte | 
Ne ae CoH aa ees 95 0 Bre 0 is adopted by the American Ephemeris (5), and appears in all editions.) 
ae a momer is concerned with practical situations and must continually endeavor to refine the fundamental constant 
which his theories rest. He therefore carries doubtful digits for the sake of consistency with other more precisely known const 
and in order to guard his numbers against end-figure error in numerical work a 
> When P, the constant of parallactie i ity i itu 
, the sti parallactic inequality in the Moon’s longiti i i 
, ality ‘ gitude, has been found by observation is case sti 
tions, the solar parallax is determined by the relation i 
ee 
P= 141566, 2 = eo 
+ me 
In Brouwer’s derivation, Hinks’ value for m@~! was adopted 
° A differential correction process based on a least squares inversion has 
astrodynamic constants by the referenced investigators 
d . ] &, 
de Vaucouleurs (30) considers me ' trig j 
purs (¢ sider: ‘an of trigonometric methods si a 5 
: : Tie ot; Ss since 1895 (8"798 | 6 : i r 
re eee al ¢ S./98,6 + 0.001,5); mean of dyn 
(8° 798,3 + 0.000,3); parallactic inequality of Moon (8"796,7 + 0.001 8); Lunar inequality of Bar os Oe 
tion (8"797,7 + 0.013,2); constant of aberration geometric (8”797 9 wi: 0 001 8) var eae ae 
ie BeOS SS MUTIE (O. 097 ,f . ; constant of aberrati spee i 
0.002,2); and mass of Earth and Moon, Pioneer V (8"797,1 + 0.000,8) of sages ee seehnei 
O.bdl, e -UUU, e 


been utilized to improve the solar parallax and o 
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TABLE 2 
Inaccuracies Due to Solar Parallax Uncertainty 


lV’ 0.0001 ; 0.001 

_ Mars 8,400 km 84,000 km 
= 5,200 mi 52,000 mi 
Venus 3,500 km 35,000 km 
; 2,200 mi 22,000 mi. 


These values when translated into correction fuel require- 
nts take on a considerable significance. 


ere m; is the Earth’s mass, the mass of the satellite 
ing negligible; r is the radius vector of the satellite’s 
mt in units of Harth’s equatorial radii, a, and 6 
declination; and the J. --- Js are the coefficients 
the higher zonal harmonics. Relations between J» , 
and J,, and J, H, and K are as follows: 


J, = 3d, J; = 3H, J, = —i5K, (7) 


ectively (Jeffreys’ D = —33J,4). 

Constants involved in the Earth’s gravitational 
ld also include 

a. the equatorial radius, 

*, the geocentric gravitational constant, 

zs ey the flattening or ellipticity of the me- 


= ridian, 

ge the acceleration due to gravity at the equator. 

« the coefficient of the term, —a, sin’2¢’ (where 
¢ is the geocentric latitude) expressing the 
depression of the spheroid or equipotential 
surface below the reference ellipsoid, reaching 
a maximum at 45°. Cf. de Sitter 1924. Its 
numerical value is variously given between 
0 and 6.8 < 10 ’. Cf. Herrick, Chapter 3 of 
Astrodynamics. 

Jacchia (47), 1958, obtained the following value 

ma study of the motions of 1957 6, and 1958 6» 


x = 3(43) x 107 
sistent with 
298.28 + 0.11 
J = 0.001624 + 0.000001 
OF EL iy | 


— 

— 

Ser) 
l| 


xX 
] 


hits of close satellites are subject to large secular 
‘turbations in the node and perigee as well as to 
ispicuous long periodic terms in other elements 
wtably the eccentricity) that are proving to be of 
at value in more exact determinations of the J; 
Equation 6. These newly determined J; are related 
E Oye q “) 
fx and @ (= a,w/k, ) as follows: 


ioeeraep— af 30 af@ + ark 
f? 4f~ 32 
f= a) yO — aK 
eliminating xk, 


f—(4+a)f+460+3h+3J,=0 (8) 


Where the values for w, the Earth’s angular velocity, 
and @ are 


2 
w 


5,317.49 X 10” radians?/sec” 
& = 0.003,461,39 


I 


corresponding to the value a, = 6,378,150 meters 
(ef. table 5). 

In Tables 3 and 4 are listed determinations of 1/f 
and Jy +--+ J¢ coefficients of the higher harmonics; 
Table 5 shows values of a,; while Table 6 presents a 
summary of the adopted geocentric constants. Note 
that seventh and higher order harmonics as well as 
longitude dependent terms have been omitted. As 
more data become available values of these astro- 
dynamic constants should be adopted as well. For a 
discussion of these constants see Michielsen (72) and 
Kaula (52). In table 6 the values of &, and g. are chosen 
to be consistent with the adopted values of a, and 1/f; 
and J. ---Js5 are consistent among themselves and 


with : (ef. Eq. 8). 


ke = 1.996,530,8( + a’ + g’/24+f') X 10° 
and 
Je = 978,036.8(1 + g’ + 3f’) 
where 


1 6.018,210.— ae. Ae 
g4= ae 
ae 2 


eet 
a 


and A = the atmospheric mass = 0.88 X 10-*. See 
Herrick, Baker, and Hilton (39) 1958. 

The J;, coefficients of the zonal harmonics in the 
Earth’s gravitational potential given in Table 4, 
have recently been derived from the observed motions 
of artificial satellites: e.g., 1958 Beta 2, 1959 Iota 
and 1958 Eta (Kozai), 1958 Beta 2 (O’Keefe, Eckels, 
Squires), 1957 Beta 1 and 1958 Beta 2 (Jacchia), 
1957 Beta, 1958 Beta 2 and 1959 Iota (King-Hele) 
etc. O’Keefe and Hekels (82) found a perturbation in 


=e = 


TABLE 38 
Flattening of the Earth 


1/f Author 


298.38 -— 0.07 O’ Keefe (79) 1958 

298.28 + 0.11 Jacchia (47) 1958 

298 .82 + 0.05 Lecar, Sorenson, Eckels (64) 1959 

298 . 24 O’Keefe, Eckels, Squires (81) 1959 
298.20 + 0.03 King-Hele, Merson (56) 1959 

298.24 + 0.02 King-Hele (55) 1960 

298.380 + 0.03 Kozai (48) 1960 (calculated from his J; 
by formula (8)) 


298.24 + 0.01 Kaula (53) 1961 

298.2 + 0.1 de Vaucouleurs (1961) (80) Weighted 
mean 

298.30 + 0.05 Adopted® 


« Adopted since it encompasses most of the more recent 
determinations and is consistent with Kozai’s values adopted 
for the J’s. 
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TABLE 4 
Coefficients of Gravitational Harmonics (all values X 107) 


ip We if Je Author 
de Sitter (25) 1924 
ee + 0.7 de Sitter (26) 1927 
1094.3 + 4.4 Jeffreys (49) 1936 
1094.08 + 0.65 de Sitter, Brouwer (27) 1938 
1091.8 + 4.2 Jeffreys (51) 1941 
1094.38 + 2.4 Spencer Jones (98) 1941 
1098.21 Clemence (19) 1948 
1089.7 == 2.1 Jeffreys (50) 1948 ~~ 
—-2.9 +0.8 O’Keefe, Eckels (82) 1958 
1082.1 —2.20 + .07 —2.15 Kozai (58) 1959 
1082.1 + 0.4 —2.359 + .002 Lecar, Sorenson, Eckels (64) 1959 
1082.53 —2.4 +0.3 | —1.7 -—0.1 +.1 O’Keefe, Eckels, Squires (81) 1959 
—2.32 Cohen, Anderle (21) 1960) € 
1082.79 + 0.15 —-1.4 +0.2 0.9 +0.8 King-Hele (55) 1960 t 
—2.6 —0.57 Michielsen (71) 1960 ¢ 
1082.190 + 0.023 | —2.29 + 0.02 | —2.13 + .04 —0.23 + .02 Kozai (to be published) 1960 (3 sat 
lites) i 
1082.28 + 0.03 | —2.29 + 0.02 | —2.12 + .05 | —0.23 4 .02 Hilton recalculation of Kozai valu 
consistent with Table 6 (43) 19 
1082.66 —2.5 —1.7 +0.3 +0.7 Michielsen (72) 1961 (perigee radi 
determination) : 
1082.61 + 0.06 —2.05 + 0.1 —1.43 + 0.06 | —0.08 + 0.11 | +0.20 + 0.05 | Kaula (53) 1961 
1082.28 + 0.3 —2.3 = 0.2 —2.12 + 0.5 —0.2 + 0.1 +1.0 + 0.8% | Adopted (1961)> 


4 J, is not necessarily consistent with the other adopted values, as they have been calculated assuming no sixth harmonic. 
b Adopted from the Hilton revision of Kozai’s values, which were based upon observations of three satellites. 


TABLE 5 
Equatorial Radius of Earth 


de (meters) Author Comments 
6,378 ,301 Clarke (15) 1880 See (86) p. 312 
6,378,298 + 34 Ledersteger (65) 1951 
6,378 , 228 Hayford Revision (93) 1953 
6,377,879 + 357 Hirose (45) 1955 
6,378,250 + 95 “Hough Ellipsoid’’ (13) 1956 Free Air 
6,378,240 + 100 “Hough Ellipsoid’’ (13) 1956 Free air with Jeffreys’ components 
6,378,285 + 100 “Hough Ellipsoid’’ (13) 1956 Isostatic 
6,378,270 + 100 Herrick, Baker, Hilton (389) 1958 i — 92970 
6,378,145 + 50 Herrick, Baker, Hilton (39) 1958 Corrected for 1/f = 298.3 + 0.1 
6,378,388 Yaplee, Bruton, Craig, Roman (110) 1958 Radar echoes from Moon 


6,378,175 + 20 


6,378, 200 Fischer (1959) 
6,378,200 + 30 
6,378,163 + 21 
6,378,255 + 35 


de Vaucouleurs (30) 1961 
Kaula (53) 1961 
Yaplee, Bruton, Miller (111) 1960 


6,378,150 + 50 Adopted (1961) 


the eccentricity of 1958 Beta 2 having a period equal 
to that of the revolution of the line of apsides, which 
theory predicts if a third order term in the potential 
with coefficient J; is assumed to exist. This indicates 
that to the third order at least, the Earth is somewhat 
“pear shaped,” i.e. the northern hemisphere has a 
longer axial semidiameter and a shorter mid-latitudinal 
semidiameter than does the southern hemisphere 
(ef. 8. W. Carey, “North-South Asymmetry of the 
Earth’s Figure,” Science 130, 1960, p. 978; O’Keefe, 


P ms 5 : ere 
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Yaplee, Bruton, Craig, Roman (110) 1958 


| Recomputed by C. G. Hilton with revised lur 
parallax and flattening 

World wide areas, recomputed by Cook (22) wi 
1/f = 298.2 

Weighted mean adopted 


Radar echoes from moon with 1/f = 298.25 + 0 
and me/meE = 81.3875 + 0.026 


Kekels and Squires (79), and also King Hele and M 
son (56)). 


3. Selenocentric Constants 


3.1. Lunar Mass. The ratio of the Moon’s mass 
that of the Earth (me) is related to the ratio of t 
solar and lunar parallaxes through the expression 


Lees He ( 
1+ me sin Te 


TABLE 6 
Adopted Geocentric Constants (1961) 


= 6,378,150 (1 + 11 X 10) 
meters 

e = 1.996,503,0 (1 + 11 X 10%) 

xX 10% megameters”? 


Adopted from Table 5 


Herrick 1958, corrected 
for new 1/f 


c second 

le = 9.780,320 (1 + 3 X 107) Herrick, Baker, Hilton 
4 meters/sec? (39) 1958 

i = 298.30 + 0.05 Adopted from Table 3 


2 = (1082.28 + 0.3) X 10-6 
3 = (—2.3 + 0.2) X 10-6 
fe (—2.12 + 0.5) X 10-7 
m= (—0.2 + 0.1) x 10-8 
6 = (1.0 + 0.8) X 10% 


Adopted from Table 4 
Adopted from Table 4 
Adopted from Table 4 
Adopted from Table 4 
Adopted from Table 4 


TABLE 8 


Distance of the Moon 
Distance (km) Author 
384,400 + 4 Christie and Gill (14) 1911; Brown (12) 


1919; Lambert (63) 1928; Am. Ephm. 
(5) 1962. 

384,407.6 + 4.7] O’Keefe and Anderson (80) 1952 

384,403 + 1 O’Keefe 1958 

384,402.0 + 1.2] Yaplee and associates (111) 1960, using 
equatorial radius of Moon = 1740 + 
1 km 

384,400.9 + 0.7} de Vaucouleurs (30) 1961 using me = 
81.36 + 0.02 dynamical parallax = 
3422"524 + 0.035 for 1/f = 300 
384,402.0 + 1.0) de Vaucouleurs weighted mean 1961 


TABLE 7 384,402 + 2 Adopted (1961) 

Determinations of the Moon’s Mass (Earth’s mass = 1) 

, Constant TABLE 9 

eee ees Reci pro cal of Moon's ae Miscellaneous Selenocentric Constants 

from Parallax 
bservation ite if Author 
87803 | 81.53 Hinks (44) 1909, Am. 0.67 + 0.03 | 1° 32’ 00” + 14” | I. V. Belkovich (8) 1949, 

Ephm. (5) 1960 on Crater Mésting A 

6.428 ,3 8.803 81.530 de Sitter, Brouwer 0.65 + 0.045) 1 32 04+15 A. A. Nefedjev, (75) 1950 

- (27) 1938 on Crater Mésting A 

6.404,9 8.803 81.829 de Sitter, Brouwer 1 33 50+ 19 C. B. Watts (105) 1955 

; (27) 1938 0.85 + 0.03 | 1 32 49+ 12 A. A. Yakovkin (109) 

6.406,8 8.803 | 81.805 de Sitter, Brouwer 1950, from Banachie- 
(27) 1938 wiecz 1910-1915 obser- 

6.439 ,0 8.790 81.271 + 0.021 | Spencer Jones (98) vations 

1941 0.60 1 3i 22 Weimer, 1954 

6.436 ,7 8.798,4|) 81.375 + 0.026 | E. Rabe (85) 1949 1 32 06 Hayn (87) 1907, p. 499, 
from Eros Am. Ephm. (5) 1962. 

6.442 ,8 8.790 | 81.222 + 0.027 | Delano (23) 1950 from 0.75 1 32 20 Hayn; cf Weimer 1954. 


right ascensions of 
Eros 
Delano (23) from 
declinations of Eros 
de Vaucouleurs (380) 
adopted 1961 
Adopted (1961)# 


5.443 ,0 8.790 | 81.219 + 0.030 


5.438 ,5 8.798 | 81.357 + 0.02 
0.001,5 
81.35 + 0.050 


4 Rabe’s value still appears to be superior. 


here L is the coefficient of the ‘lunar equation,” 
hich is involved in the reduction of geocentric co- 
dinates to the center of gravity of the Earth-Moon 
stem (barycenter). Several determinations of the 
oon’s mass (in units of the Earth’s mass) were based 
the close approach of 433 Eros in 1930-31. Assuming 
e solar and lunar parallaxes to have been found 
eviously, the method consists in comparing the 
served positions of Eros when nearest the Earth 
th an accurate ephemeris, fitting the residuals to a 
00th curve that has the periodicity and zero points 
the lunar equation, and using the curve to improve 
e adopted value of LZ. Some results are shown in 
ible 7. Improved values of the Moon’s mass will no 
ubt be obtained from lunar satellites. 

22. Lunar Distance. Until recently the accepted 
lues for the Moon’s distance and parallax stemmed 


from meridian observations of the Moon’s declination 
made during 1906-1910 from Cape of Good Hope 
and Greenwich Observatories by Christie and Gill 
(14). The reduction of their measures, which required 
an assumption as to the ellipticity of the meridian 
between the two stations, gave the following results: 


Moon’s parallax = 57’ 02”70 
“constant of the sine parallax” 
= 3,422”54 (ie. sine p/sine 1”) 
distance = 384,400 km. 


These values are still published in the American 
Ephemeris (1962, p. 485). 

Lambert (63) in 1928 reconsidered the reduction 
of the observations by deriving geometrical lunar 
parallaxes for values of the ellipticity of the meridian, 
1/f, ranging from 292 to 300, and compared them with 
equivalent values of the dynamical parallaxes, using 
as constants® 

6 This value of g, though somewhat higher than accepted 
values, was considered appropriate in view of gravity meas- 
urements from a submarine made by Vening Meinesz (1927) on 
a voyage from Holland to Java, which indicated that gravity 
at sea slightly exceeds gravity on land. 
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a, = 6,378,388 km (International Ellipsoid) 
978,052 cm/sec” 


© 
Il 


I 


yw = 81.53 (Hinks’ value from Eros). 


Lambert found that for geometrical and dynamical 
parallax equal to 57’ 02”72, 1/f = 293.5. For the more 
acceptable flattening of the International Ellipsoid 


(297) a discrepancy of 0’16 might be interpreted 
due to inexact constants used in computing the d 
namical parallax, or to errors in the deflection of f 
vertical, the elevation of the geoid above the sphera 
or the observations of Crater Mésting A used in ¢ 
riving the geometrical parallax. He therefore retain 
Christie and Gill’s value until more observational d 
from other meridians should become available. T 


TABLE 10 
Reciprocal Masses of the Planets (Sun’s Mass = 1) ; 
Planet mo Author Method : 
Fl 
Mercury 6,000,000 Newcomb (76) 1895 v 
(Adopted by Am. Ephm. (5)) 2 
Clemence (16) 1953 i 
7,500,000 + 1,500,000 de Sitter, Brouwer (27) 1938 weighted mean 4 
6,120,000 + 43,000 Rabe (85) 1949 Eros 5s 
5,970,000 -+ 460,000 Duncombe (81a) 1956 Venus - 
6,100,000 + 50,000 Adopted (1961) A 
3 
Venus 408 , 000 Newcomb (76) 1895 ; 
(Adopted by Am. Ephm. (5)) ; 
406,358 + 723 Fotheringham (34) 1926 Earth, Mars and Mercury 
403,490 -+- 2,400 Ross (90) 1916 Mars 
404,700 + 800 Spencer Jones (97) 1926 Sun 
404,000 + 1,000 de Sitter (27) 1938 weighted mean 
407,000 + 500 Morgan, Scott (73) 1939 Sun 
409,300 + 1,400 Clemence (18) 1943 Mercury 
406,645 + 208 Rabe (85) 1949 Eros 
407,000 + 1,000 Adopted (1961) 
Earth + Moon 329,390 Newcomb (76) 1895 
(Adopted by Am. Ephm. (5)) 
327,900 + 200 de Sitter (27) 1938 | weighted mean 
328,390 + 108 Witt (106) 1933 | Eros 
328,452 + 43 Rabe (85) 1949 | Eros 
328,446 + 48 E. Rabe (85a) 1954 Eros revised for precession 
328,440 + 40 de Vaucouleurs (30) 1961 | Adopted 
328,450 + 50 Adopted (1961) | 
Mars 3,648 ,000 Leveau (67) 1890 Vesta 
3,098 , 500 Newcomb (76) 1895 
(Adopted by Am. Ephm. (5)) 
3,601,280 Leveau (66) 1907 Vesta 
3,085,000 + 5,000 de Sitter (27) 1938 weighted mean 
3,110,000 + 7,700 Rabe (85) 1949 Eros 
3,079,000 + 5,702 Urey (1038) 1952 Deimos 
3,090,000 + 10,000 Adopted (1961) 
Jupiter 1,050.36 Eneke (83) 1938 Vesta 
1051.42 Hansen (86) 1865 Kgeria 
1,047 .538 Kriiger (60) 1865 Themis 
1,645.63 Leveau (67) 1890 Vesta 
1,047.855 + 0.065 Newcomb (76) 1895 weighted mean 
(Adopted by Am. Ephm. (5)) 
1,047.34 Newcomb (77) 1895 Polyhymnia 
1,046.04. Leveau (66) 1904: Vesta 
1,047.40 + 0.03 de Sitter (27) 1938 weighted mean 
1,047.4 + 0.1 Adopted (1961) 
Saturn suey ' Newcomb (76) 1895 weighted mean 
3,490 + de Sitter (27) 193 1 
3,497.64 + 0.27 Hertz a Hig ihe ge sc 
ee 0.4 eae (16) 1953 and (17) 1960 Jupiter perturbations 
) dopted (1961) 
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TABLE 10—Continwed 


Planet mol 


Author Method 


22,869 


22,750 + 200 
22,800 + 100 


_ Neptune 


19,314 
19,700 
19,500 + 200 
19,500 + 200 
Pluto 332,488 + 76,472 


360,000 
307 ,000 
350,000 + 50,000 


St 


Newcomb (76) 1895 
(Adopted by Am. Ephm. (5)) 
Clemence (16) 1953 

de Sitter (27) 1938 

Adopted (1961) 


Clemence (16) 1953 

(Adopted by Am. Ephm. (5)) 
Newcomb (76) 1895 

de Sitter (27) 1938 

Adopted (1961) 


Nicholson and Mayall (78) 1931 
Clemence (16) 1953 

Kuiper (62) 1950 

Adopted (1961) 


@ Actually given as 1.08 X Earth’s mass + 0.23. In Pub. A.S.P., 43, no. 74, 1931, the authors give Pluto’s mass as 2 the Harth’s 


ass. 


compares the older results for the Moon’s distance 
ith those of O’Keefe and Anderson (80) derived from 
n0toelectrically recorded star occultations (1952) 
id with the results of radar reflection observations 
y Yaplee, Bruton, Craig, Roman and Miller (1957— 
160). 
3.3. Figure of Moon. The Moon’s figure appears to 
» best approximated by a triaxial ellipsoid with 
miaxes a, directed toward the center of the Earth, 
coincident with the rotational axis, and 0b, per- 
ndicular to a and c. Alexandrov’s (4) values derived 
om Yakovkin’s data on the Moon’s physical libra- 
mn, given below, indicate that it is immaterial whether 
e librations are assumed to be forced or free: 


forced libration free libration adopted (1961) 


semiaxis a |1738.67 + 0.07/1738.57 + 0.07/1738.57 + 0.07 
semiaxis b {1738.21 + 0.07/1738.31 + 0.07|1738.21 + 0.07 
semiaxis c |1737.58 + 0.07|1738.58 + 0.07|1737.58 + 0.07 


loments of inertia about the respective axes: 


mee (hy? +c’), B= © ke +a), 
5 5 


(10) 


C =e (a +0’) 


sre derived by Alexandrov as follows: 


forced libration free libration 


A X 10% gm/km? 88.837 + 0.024 88.838 + 0.024 
BX 10° gm/km? 88.856 + 0.024 88.856 + 0.024 
CX 10% gm/km? 88.893 + 0.024 88.893 + 0.024 


3 


Jeffreys (48) finds the following relations: 


C—A 
goes 0.000,626,6 + 0.000,002,7 (standard error) 


dle Je 
C 


= 0.000,204,9 + 0.000,000,9 for a forced libration 


= 0.000,209,8 + 0.000,002,2 for a free libration. 


3.4. Miscellaneous Selenocentric Constants. If a fune- 
tion of the moments of inertia about the principal 


axes is defined by 
CB 7B 
yn (stam = Po 
foal ees a) 


its value and that of the inclination of the Moon’s 
equator to the ecliptic can be found from observations 
of the physical bration. Following in Table 9 are some 
recent determinations from the Transactions of the 
L.A.U., 1952 and 1955. 

H. Hirose and R. Manabe (45), from observations 
of occultations of the Pleiades, applying Hayn’s cor- 
rections for irregularities of the limb, obtained a value 
for the Moon’s semidiameter of 


039780. 0.07, 


The value of the ratio of the Moon’s diameter to that, 
of the Earth, used by H. M. Nautical Almanac office 
in occultation predictions, 1.e., 


k = 0.272,495,3 
gives a semidiameter of 1,738,016 + 14 meters, as- 


suming a, = 6,378,150 + 50 meters. 


4. Planetocentric Constants 


Heliocentric constants are important over a major 
part of an interplanetary trajectory, but the terminal 
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phase is influenced by such constants as the mass, 
diameter, and flattening of the target planet. Plane- 
tary distances in astronomical units and geocentric 
and heliocentric coordinates are generally well es- 
tablished. Planet masses enter into the equations of 
motion of a vehicle through the perturbative terms. 
Older values of the masses, found from planetary 
satellites, the perturbative action of one planet on 


another, or on a minor planet or satellite, are in the | 
cess of refinement to meet the need for greater 
curacy. In his definitive differential correction of 
orbit of Eros, Rabe (85), for example, derived | 
provements to the masses of the four inner planets 
well as to precessional and other constants and cert 
elements of the Earth’s orbit. Similar improveme 
of planetary masses will be derived from observatic 


TABLE 11 = 
Planetary Diameters (Angular Diameters Are Reduced to 1 a.u. Unless Otherwise Stated) ~~ 
Planet Linear Diameter (km) Angular Diameter Author 
Mercury 4,842 6”68 Leverrier (5) 1843 
(Adopted by Am. Ephm.) 
4,670 6.45 Dollfus (88) 1953 
4,650 6.42 Muller (88) 1953 
4,676 + 31 de Vaucouleurs (31) 1961 
4,660 + 30 Adopted (1961) 
Venus? 12,191 16.82 Auwers (5) 1894 
(Adopted by Am. Ephm.) 
12,640 W. Rabe (87) 1928 
12,513 Ross (91) 1928 
12,060 16.64 Muller (74) 1948 
12,246 Kuiper (61) 1949 
12,200 + 10> Menzel, de Vaucouleurs (70) 1960 
12,200 + 20 Adopted (1961) 
Mars 6,743 + 22 9.268 See (94) 1901 
9.222 (polar) 
6,784 9.36 Hartwig (5) 1879 
(Adopted by Am. Ephm.) 
6.652 + 11 9.175 + 0.015 Trumpler (101) 1927 
6,577 + 34 (polar) 
6,679 + 42.3 9.48 + 0.006 Van de Kamp (104) 1928 | 
6,860 + 21.7 W. Rabe (87) 1929 . 
6,834 + 50.2 9.44 + 0.07 Reuyl (89) 1941 
6,826 9.29 (polar) Muller (74) 1948 
6,830 de Vaucouleurs (31) 
6,830 + 10 Adopted (1961) 
Jupiter 142,745 196794 Sampson (5) 1910 
133,286 (polar) 183.82 (Adopted by Am. Ephm.) 
142,750 + 100 Adopted (1961) 
Saturn 120,798 166.66 H. Struve (5) 1898 
108 ,096 (polar) 149.14 (Adopted by Am. Ephm.) 
121,000 + 100 Adopted (1961) 
Uranus 49 ,693 68.56 weighted mean, Barnard (5) 1896, See (5) 1§ 
Wirtz (5) 1912, 
(Adopted by Am. Ephm.) 
49,700 + 100 Adopted (1961) 
Neptune 52,999 eave Barnard (5) 1902 
' (Adopted by Am. Ephm.) 
ees 600 2.044 at 30.07 a.u. Kuiper (62) 1949 
50,000 + 500 Adopted (1961) 
Pluto 0.46 as diameter (0.23 + 0.01 at 35.56 a.u. Kuiper (62) 1950 
a + 1000 Adopted (1961) 


* Diameter of atmospheric cloud level shell. 
» Angular diameter up to half-intensity point of 
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atmosphere on measurement of occultati 
and linear diameter 12,321 + 8 km. They estimate this half-intensity cm, above Rede 


the linear diameter of the cloud shell is approximately 12,200 + 10 km 


point to be 65 km. above the cloud level (+ 10 km) so t 


. ring close approaches of Icarus (88), Griqua (86), 
iberga (11) and several other minor planets which 
ave been discovered in recent years. 

| Table 10 contains the principal determinations of 
lanetary masses, including the atmospheres and 
atellites where these are present. The diversity in 
he values of linear planetary diameters given in Table 
1 illustrates the difficulty of obtaining precise measures 
f the angular diameters (reduced to a distance of 1 
u.) especially in the case of planets with atmospheres. 
lore accurate determinations are greatly needed for 
adar measurements of distance and, particularly, 
hen landings are to be attempted, since positional 
ata always refer to a planet’s center. 

The flattening of a planet, defined as (a — b)/a 
rhere a and b are the equatorial and polar radii re- 
pectively, and the gravitational harmonies, /;, be- 
ome important when orbital motions of nearby satel- 
tes are under consideration. They may be found 
rom dynamical relations, assuming hydrostatic equi- 
brium and a particular density distribution, from 
he motion of a planetary satellite, while the flattening 
ay be measured on photographs directly. Table 12 
ontains some estimates of the flattening of Mars, 
able 13a values for the flattening and J; of Jupiter, 


TABLE 12 
Flattening of Mars 
1/f Author Method 
Struve (99) 1895 dynamical 
Woolard (108) 1944 dynamical 
Trumpler (101) 1927 photographs 
de Vaucouleurs (29) weighted mean of opti- 
1950 cal measures covering 
100 years. 
77 + 30 | Dollfus (88) 1953 optical in yellow light 
Muller (88) 1953 optical 
Camichel (88) 1954 optical 
0 + 50 Adopted (1961) 


a Adopted in order to encompass more reliable dynamical 
eterminations and Trumpler’s photographic determination. 
‘Keefe and Kopal have suggested in conversation that the 
iscrepancy between optical and dynamical determinations 
ay be due to the depression of the tropopause in the Northern 
egions of Mars. 


TABLE 13a 


Flattening and Gravitational Harmonics of Jupiter 


1/f J» J4 Author and Method 
15.0 Sampson (5) 191C, 
: from polar and 
| equatorial radius 
15.2 0.014,93 |—0.000,405 | Sampson (92) 1921, 


from satellites 
0.014,90 |—0.001,0 Brouwer (10a) 1946 
Deo 0.014,71 |—0.000,336 | De Marcus (24) 1959 
0.014,8 —0.000,35 | Adopted (1961) 
+0.000,1} -+0.000,05 


TABLE 13b 
Flattening and Gravitational Harmonies of Saturn 
1/f£ Je Js Author 
10.2 1.016,67 | —0.000,515 | DeMarecus (24) 1959 


0.016,107 | —0,000,865 
+0.000,007 | +0.000,026 | Kozai (59a) 1957 


10.2 0.016,10 | —0.000,865 | Adopted 1961 
+0.000,01 | +0.000,026 
TABLE 13¢ 


Adopted Values of Gravitational Harmonics and Flattening 
of the Planets 


Planet 1/f Jo Js Comments 


Mars 150+50 | 0.002,011 1/f from Table 
12 and J» from 


De Marcus 


(24) 1959 
Jupiter | 15.2+0.1) 0.614,8 |—0.000,35 | Adopted from 
=+0.000,1 |+0.000,05 Table 13a 


Saturn 10.2 0.016,10 |—0.000,865) Adopted from 


+0.000,01/--0.000,026) Table 13b 
Neptune | 58.5 0.004,9 De Marcus (24) 
1959 


Table 13b values for the flattening and J; of Saturn, 
and Table 13¢ some adopted gravitational constants 
for four planets. 


o. Atmospheric Constants 


Experience with geocentric artificial satellites has 
already demonstrated the importance of atmospheric 
constants to astrodynamics. In fact the atmospheric 
densities at high altitudes have required revision 
of an order of magnitude over older accepted values in 
order to agree with satellite drag perturbations. Because 
the Earth’s atmosphere has been discussed in detail 
in various sources (cf (6) and (7) for example), the 
following section will treat primarily the atmospheres 
of Venus and Mars as being of special immediate in- 
terest to the space traveler. Atmospheric constants 
include composition and density at surface and their 
variation with height; diurnal, seasonal, solar and 
latitudinal fluctuation in density; and winds. As in- 
terplanetary probes make available improved as- 
tronomical data, both the quantity and accuracy of 
atmospheric characteristics will undoubtedly improve. 

The variation of density with altitude is best defined 
in terms of equation (57) from (7), Le. 


: * 
inp = inp, — (1 ei 


Lin x \ 
{in 1 + (T.,.)p (H A, | 


where the geopotential altitude for a spherical planet is 


= Ago 
Go (+> H/ i’) 


H* (13) 
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atmosphere is subject to forces and stresses at 
altitudes that result in complicated flow patter 
Hence rather large scale unpredictable variatior 
from any model would be the rule, not the exception 
At this time, if an adoption of a wind model is to } 
made, the summary work published by J. M. Sok 
(95) is recommended. 

The seasonal, diurnal, solar, latitudinal, ete. va 


Thus, given values for the base altitudes, H,* ; the 
base molecular scale temperature, (7'm)»; the base 
mean molecular weight, m,; and the temperature 
eradient, L,, ; from Table 14, an atmospheric density 
profile can be constructed. The sea level composition of 
planetary atmospheres is given in Table 15. 

The question of a model atmospheric wind pattern 
is almost meaningless even for the Earth since the 


TABLE 14 


Atmospheric Structure Model 


Temperature gradient 
Planet Hy* (km) (Deck) mb P;* (Newton/m?) po (kg/m) Ln i jo ae 
Venus™ 0 340 + 15 42 .4¢ (2.2 +1) X 10% 4+3 —5 0<H<20° FF 
(580) (10 to 30 X 105)* : 
(535-675)! 
(350)¢ 
go* = 8.42 m/sec? | 18 + 3 250 + 30° 42.40 (1.7 + 2) X 10% ey iw) +10 20<H<30° — 
(cloud level)| (285) * Gbeio6 N03) -+-9k i 
R’ = 6,100,000m | 70+ 10 | 240 + 204] ~40 ? 2 1.5 70<H<804 — 
85 + 10 297 + 10° ~A40 0.26 + 0.18° (4 + 1) X 10° +3 80<H<100° 
Earth: 0 288.16 | 28.966 1.01 X 105 1.2250 6.5 O0<H<ll | 
go* = 9.80665 ila 216.66 28 .966 228 >< 104 3.6480 X 1071 0 11<H<25 : 
m/sec? : 
R’ = 6,356,766 m 25 216.66 28 .966 2552) < LOs 4.0639 X 107 =+- 3.0! 25°< Hm 4 
47 282.66 28 .966 1.25 X 10? i aayin s< Ur 0 4] <= ASS 
53 282 .66 28 .966 6.14 K 161 (eT NU —4.5 53 < Hs 
79 165.66 28.97 iL 2B} 202A) eat Ome 0 19< H< 
96 165.66 28.97 ieee) >a 2.846 xX 10°& +4.0 90< H< 
105 225.66 28.86 9.66 X 107? 1-542 x 10s +20 105< H< 
160 a20.06 28 .04 A 10 Xo =* Len 23 Ee Om +10 l160<H< 
170 1425.66 Pit NG) SLA Oe (shy. Seis +5 170<H< 
200. 1. 575.66 26.328 1.65 xX 10° Saye! Ce +3.5 200 < H < 70 
Mars 0 290 + 15° 28 3° (8.7 + 0.4) X 10% | (0.99 + 6.05) X 107 —-4 0Q0<H< 
(218 + 76°)!| (44 to 92) 
Got = 3.91 m/sec” 30 + 10° 180 + 20° ~28 (1.4 + 0.2) X 108 |-@.7 £1) X 107 —=1 “30 <—Hea ac 
R’ = 3,400,000 m 90 + 205 140 + 10° ~28 7.2+2 (E> Sea ee dOiet +1.5 90 < H< 110 
(130)" oon 
200 + 50! 300 + 50! ~28 ~!1 ~1 & 107-5 0 ae 
350 + 50) | 300 + 50) | 288tol4i | ~1 x 107 ~1 X 107 ~50 350 < H < Sap 


. a ea = | bar. The surface pressure is determined in the case of Mars by molecular light scattering. De Vaueor 
eurs 960 finds Pr = (8.5 + 0.4) & 10° nt/m?, while Urey (102) 1959 modifies Dollf 1957 : eo 
) ; By be iS} s? 19 § 08 
b pb = Pymp/ (Tm) R*. : < - te me is a : ts hh 
¢ Urey (102) 1959; approximate values for the uncertainties added 
| ; xime alues rte 3g d by the authors; pressure of 1.7 C 2 < 
Venus’ surface was taken as basis for computation of ‘sea level’ values. ha amarante 
4 Pettit and Nicholson (83) 1955. 
e Differential refraction studies of photoelectric observations of i i 
‘ elect, servations of the occultation of Regulus by Venus, Me f 
(70), under the assumptions that g = 8.6 m/sec, 90% COn, 9% No, and 1% A+ --- . We assume 65 a a Mee 
7 ass : KK a Ss fe u 3 


above surface. 


f ARDC Model Atmosphere (6) 1959. 


ni 1 a 
* From an altitude of 180 km on, the molecular weight is obtained from m = 


4 Kopal (57) 1960. 


27.106 — 7.935,697,10 tan! 


(H — 180) 
140 


AIeSS Karrer INC 3S fe 5 t p 

I eSs st 1 Kiess 54 1960. ») ‘h \ molec y V =e} é { \ e e ense atmos Lere tha 
? c ( ) Suc ral € ‘ular W eight V ould Eek ult in a two to four tim mor 1 } 

has been quoted here. ’ ; fe : . 


i Hstimated from tables given by Yanow (112) 1961 for po = 107 kg/m 


HLO. 


( ) € 5 »f radio ast Ono 9, ol » ‘ t V + TI + AA Ord X or Roy ¢ 
€ ‘ ‘ J 318 © MUI IDS V 5 S ng e 2. 2, € t 3 
rape M1¢ al servations oO 0 9) nN 

k B rett (a 1961 ound on t i € pas f{ nus, assum 1g (0 /o CO, 22— q 2a d 0 ne) 


' Mayer McCullough, and Sloanaker (68a) 195 1 i 
j ; et: , ane Ot 958, find a surface temperature for r 6 = erore e 
A : : 7 , ts at > for Venus of 620 110°K 5 i 1 
i > : ie eee ; ; an ier ce : “ X 40 da "s bef 
c( njuncwon, wnd oot) + 73 ‘K at conjunction. In 1960 (68b) they found radio surface tem yeratures of 5385 + Ae < 16 >: 
conjunction and 675 + 85°K, 39 days after conjunction. : aa 2 ae dav — 


me A recent yaper by Cr yal Ss a y Le 7 i/ q 3 849. 858 we yi > t , + t a V q 1 
Pe mosp e : . ze rs Le WD. (Sc NCE 3 > oe : fe ) 21 ves an @xceé llen J summary of our presen ) knowledge of enus including 1 
2 a S] lere. nae 1 Seems to favor a hot d Dy surface with an at mospheric pressure of about 4 newt ons/m and a aera re of abot 
assuming ¢ 7 a re a 1 ‘ ; ‘t Re . 
) ga emper it ure gI adient of = 10 (0 <4 H < Be) km) and a ropospheric heigh of a yout 23 km abo ve the pl anet Ss 8 f C 
() K SS t 7 c 0 c 7 2 é »8 Surta 


ae 
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TABLE 15 
Atmospheric Composition at Surface (% by Volume) 


& above cloud layer. 
> Kuiper (61) 1949. 
© Urey (102) 1959. 


ars is composed of Nitrogen oxides; principally N»O and N2Oy. 


ion in the atmospheric densities, like the winds, are 
bt only quite variable but as yet not completely 
Mfinitized. Consequently we can only recommend to 
e reader the research in the field by L. G. Jacchia 
6) and by H. A. Martin and W. Priester (68) and 
opose no model for adoption at this time. 
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Rigid Body Attitude Stability and Natural 


Frequencies in a Circular Orbit 


D. B. DeBra and R. H. Delp’ 


Introduction 


The equations of attitude motion of a small rigid 
body in orbit have been recorded frequently in the 
literature; e.g., [1], [2]. Subject to several restrictions 
and assumptions they can be reduced to a set of three 
simultaneous linear, ordinary differential equations 
with constant coefficients. The conditions are that: 
(1) the only forces on the body are representable by 
an attractive inverse-square central-force field about 
an unaccelerated point; (2) the body is small enough 
so that attitude motions have no significant effect on 
the orbital motion of the center of mass; (3) the orbit 
is circular; (4) the attitude deviations from the equilib- 
rium position are small. The equations given in [2] 
are based on these assumptions. They describe the 
attitude behavior of a rigid body subject to gravity 
gradient torques in a circular orbit of angular speed wo . 


Stability 


With the coordinates of Figure 1 the equilibrium 
positions correspond to the body principal axes aligned 
with the coordinates. The principal moments of inertia 
about these axes are [; , J, , /; respectively. 

From [2] the equations for small motions about this 
equilibrium position are: 


1,6. + (Is — In)ao 01 — wb + In — Is) = 0 


wo (Ty + I, — Is) + Tob. + (13 — I) Op Go 


») 1) 
= —30) (J; — lh) 


9 


1:63 = —3 wo (Lo — I,) 63. 


Taking the Laplace transform; replacing the frequency 

. 2 e 4 nea 
variable s by p = (s/wo); introducing moment of 
inertia parameters by setting: 


(Iz — I2)/I, = ky, (I) — I) [fo = ke 


and (I, — LEQ a i =k; 


bo 
= 


. "poe 3 . . . A 
and simplifying, the equation can be written in matrix 
form as: 


1 Received October 1960. 

? Lockheed Aircraft Corporation, Sunnyvale, California. 

* When this study was made in 1958, the moment of inertia 
ratios I,/I; and I2/I,; were used [3]. In March 1960 R. Smelt 
introduced the notation of (2) in his course “Dynamics of 
Space Vehicles” at Stanford University. Smelt’s notation con- 
siderably simplifies the analysis. 
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Bae -Vpi-h) 0 | 


Vp (1 + ke) (p — 4k) 0 
if 0 0 (p + 3k) | 
6; (p) | 0 
05 (p) =| 0}. 
63 (p) 0 


Since the equations are homogenous, the behavior 
obtained from the characteristic equation which 1s jus 
the determinant of the square matrix in (3): 


[p° + (1 — 38ke — kiko) p — 4kyke] (p + 3k3) = 0 


This equation is a cubic in p = (s/w). For th 
equilibrium to be stable it is necessary that all the root 
of (4) be real and negative. The conditions this place 
on the coefficients of (4) are therefore conditions on th 
shape of the body for stability to exist. Together wit 
the physical limitations on the k’s, a chart can b 
prepared describing the necessary shape of a body fe 
stability to exist (Figure 2). These conditions are: 


Condition I 


For a principal moment of inertia where x2 and ; 
are the coordinates of an element of mass dm, I, = 


RADIUS FROM 
FORCE CENTER 


Fra. 1. Orthogonal set of unit vectors defining the local o 
bit coordinates and the small angles which define body attituc 
relative to these coordinates. 


2 2 
rate [tet = a8) am 
[ki | = i \= oe ll 


e 


f. + x3) dm. Then 
e- | (29" + 23°) dm 


and similarly —1 < ky < 1, and —1 < ks < 1 for any 
physical body. 


Condition IT 
From (4) the factored root requires k; > 0. 


Condition IIT 


Again from (4) the quadratic formula must give 
two negative real roots for p; hence: 


(A) (1 — 8k. — hike) > 0 

(B) (1 —3ho — ks)? + 16 hike > 0 
and 

(C) kak, < 0. 


Because any two of the k’s describe the shape of the 
body, one of the k’s may be eliminated. From (2) 


ky + ko + ks + ky ko kz = 0 80 that: 
ky = —(ko + ky)/(1 + keoks) 


ky = —(ks + ky)/(1 + keshky) (5) 


and 
ks = —(ky + kp) /(1 a kik). 


In view of the complexity of Condition III, k; will 
be eliminated. Since it follows from Condition I that 
1 + kike > 0, Condition II becomes kz + ky < 0. The 
stable regions in a kik, plot are therefore bounded by 
parts of the following curves: 

Condition I 


ky = —1 and k= —1, 


Condition IT 
ky = —he, 
Condition IIT 
(B) hm = — [7 = 3) 4 4n/e Dae 
(C) kp =O and) km =0. 


Condition III (A) and parts of each of the other condi- 
tions are weaker and do not contribute to the boundary. 


k= (Is-T, )/1, 


so -04 
0.8 es 
> - 
06 2 > 
0.4 | -0.5 
0.2 3 = 
Oo. 3 a 
-|.0 -O8 -06 -04 -0.2 0 x 0.6 
ky 
‘NATURAL FREQUENCY ON CURVED BOUNDARY =07 


-0.8 


Wi 


Py est 
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Fia. 2. Attitude stability and natural frequencies of a rigid body in a circular orbit plotted versus hk; and ko. 
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Fra. 3. Attitude stability and natural frequencies of a rigid body in a circular orbit plotted versus r; and re. 


The resultant regions of stability are shown in Figure 2 
and Figure 3. 


Natural Frequencies 


For each point in the stable region there are three 
natural frequencies associated with the natural motions 
of the body attitude. The motion about the 1 and 2 axes 
is coupled, so both the natural frequencies from the 
quadratic part of (4) are present in the 6; and 4 mo- 
tions. These frequencies are called w, and a». Since 
they may be incommensurable, the resultant Lissajous 
figure swept out by the tip of 3 may not be closed; 1.e., 
it can be a conditionally periodic motion. The @; motion 
occurs about the 3 axis at frequency ws . 

The natural frequencies can be obtained from (4) 
and expressed in terms of any two of the k’s using (5). 
iiep, = = (oi oule 2 = —(e/m)” and 
Ps = —(s/wo)’ and k; and kp are given: 


(io) == (-—1 -- ake +- kike 
+ -V/1( — 38ke — kik)? + 16k) /2 


6 
and 7 


D3 = 3(ky oo ke) /(1 + Kyke). 


The explicit expressions in terms of the other two com- 
binations of the k’s can be obtained by direct substitu- 
tion using (5) again. As an example Table I has been 
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prepared for two sets of values of k; and ky. The corre 
sponding points are shown in Figure 2 as point @ 
and Qe. 

Though the analysis is simpler using the / parameters 
the ratios r, = )/I3; and re = I2/I3 give a simpleé 
visualization of the physical body. It is therefore usefu 
to have the relations between the / parameters an¢ 
the r parameters. From (2), (7) gives the moment 0 
inertia ratios as functions of the k’s taken in pairs 
For example the r; and rz corresponding to points @ 
and Qs»: can be computed from (7) and are included i 
Table I. 


ik ee i + ke Bi, | ea kg x ] + ko kg 
: 1 + ky kp il +. ky 1— ko 
(7 
Se zi ‘ae Key ey il 25 ky ke te it! + kg 
é 1 4 ky ko 1 + ky 1— ke 
The inverse functions are: 
ky => (1 — ro) [Ty 
ke = (ry — | )/1e ; 
(8 


and 
ks SS 79 te 
In order to avoid resonances or to achieve certai 


dynamic behavior, constraints may be placed on tl 
natural frequencies. It is then desired to determit 


: r TABLE I 


= . GQ) (3) @ 
A ko Point kg o a as rat re 
a 3 wo ao ao 


oo-5j}—0.8 | Qi | 0.5 | 0.695 | 1.82 | 1.225 | 0.333) 0.833 
=—0.5 | 0.04 | Q» | C.469} 6.3165) 0.928) 1.187 | 1.06 | 1.53 


i 


whether these frequencies are realizable and to find the 
jody shape as described by the k’s or the r’s. 

Only two of the natural frequencies are independent, 
however, since two parameters describe the shape of 
the body. Because of the coupling of the motion between 
the 1 and 2 axes w, and w cannot be distinguished from 
ach other. There are therefore only two ways the 
requencies may be given, viz. w, and w., or w; and 
aither w, or wo. Then from (4) if o; and we are given 


feain py = — (e/a) etc.). 
ky = = 3D4 po/\pr Po + AC] = Di po) | 
kee = [pi D2 +401 + rr + pn) /12 (9) 


a oOp1 P2 [pr D2 + AQ + pi po) | 
3(4 — pi po)[pi po + 40. + pi + p»)] 


wT given w; and ws 


ky = (A + V/A? — B)/6 (p; + 4) 


ks 


vhere 
A= p(pr +-2pi + 4) — 991 
B= 36p.(m + 4)(1 + mi — Bs) 
ky = OC £.4/C? — D (10) 
vhere 


% = (4p; — 9p, — 2 pips — pr ps) /6(~1 + 3pips — 4) 
= i(p: + 1)/(pi + 3p1p3 — 4) 
ks = p3/3. 


An example is worked in Table II for a given a, and 
9; . Two possible configurations are obtained from (10) 
ind are shown on Figure 2 as points Q; and Q, . 

For even slight eccentricity there is parametric and 
orced excitation at wy. The points Q; and Q, and any 
ther point where a natural frequency is equal or near 
bit frequency are excited near resonance and would 
9 particularly poor choices. 

Under some circumstances the body configuration 
nay be a secondary consideration to the frequencies. 
tiven two natural frequencies it may be desired to 
‘btain the third frequency directly rather than by 
uccessively solving (9) or (10) and then (6). The 
xpression for the third frequency given ; and ; is: 


p, = (—2B + »/4B? — 16AC)/A (11) 
vhere 
A = pr(1 + ps) + 4p: (2 + ps) + 16 
B= pr (2 + ps) + pi + 4(2 — Ds) 
C= py + pi(2 — ps) + (1 — da) 


TABLE IT 
2 25 | Point ky ke kh sae r re 
wo Wp wo 
0.707) 1.0 | Q; 0.527 |—0.733 | 0.333 | 1.76 | 0.433 |0.768 


Q: |—0.407 | 0.0856} 0.333 | 0.524) 1.124 |1.458 


and given w; and a, it is 


[pr po + 41 Sees pe) — 36p1 pe 
(4 — pr po) [pr po + 4. + pi + pr»)] 


Though the preceding formulae are all useful for 
calculating precise results, they are slow for preliminary 
design studies. The loci of natural frequencies have 
therefore been plotted on the stable regions of Figure 2 
and 3. (These loci are for the normalized frequencies 
w/w etc. not the p’s.) Any pair of parameters (except 
w; and k;) may be used to enter the curves and locate a 
point. For example given w; and r, the designer can 
immediately read off the associated w,, w. and 7 from 
Figure 3. 


(12) 


ps = 


Conclusion 


In order that a small rigid body in a cireular orbit 
may have an orientation which is a position of stable 
equilibrium in all three axes, it is necessary that it 
have a mass distribution which is one of two types. The 
first type corresponds to the stable region in the fourth 
quadrant of the k,k. plot of Figure 2, and is the type 
usually referred to in the literature. It is characterized 
by 4, < I; < I;. The second type corresponds to a 
stable region in the second quadrant of the k, /» plot. 
In this second region the vehicles are characterized by 
Iz; < I, < I,. Vehicles of the first type have a mass 
distribution similar to a board flying edgewise along 
the orbit with its long axis vertical. The second type is 
like a board flying broadside along the orbit with the 
long axis normal to the orbit plane. The second type 
is stable because the destabilizing gravity gradient 
effect about the 2 axis is weaker than the gyroscopic 
stabilization due to angular velocity about the 3 axis. 

For the assumptions made at the beginning of this 
note not more than a total of two shape parameters 
and/or frequencies may be chosen independently. 
Any pair (except w; and k;) completely determine the 
other shape parameters and frequencies. 
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Analysis of Error Progression in Termina 
Guidance for Lunar Landing 


; 


P. J. deFries’ 


Abstract 


A lunar descent scheme with n periods of engine ON-OFF 
is investigated. Each period may have a different thrust 
level but it is constant for the length of the period. The 
engines are ignited by altitude and cut off by velocity signals. 
Five groups of errors are considered: ignition altitude (As), 
velocity change (Au), thrust level (A/’), lunar gravitation 
(Ag), and initial approach velocity (Aw). Equations are 
derived that allow the computation of the influence of any 
given error occurring during any one of the m periods on the 
final velocity at touchdown. 

It is concluded that it is principally advantageous to break 
up the descent into ON-OFF periods. Also it is inferred, as 
far as guidance is concerned, that the problem of descent 
rests with the very last period of braking. 

The upper periods can always be handled with constant 
thrust and rather crude instrumentation whereas the in- 
strumentation for the last period depends upon the required 
softness of the landing which also determines whether or 
not variable thrust is necessary. 


LIST OF SYMBOLS 


M Mass 

M Mass per unit time expelled by rocket motor 

Mo Initial mass 

i Time 

g Gravitational acceleration 

go Gravitational acceleration of moon at surface 

a Gravitational acceleration of moon at end of braking 
period 

y Momentary displacement of vehicle (vertical) 

y Momentary velocity of vehicle (vertical) 

PF Thrust of rocket motor 

In Natural logarithm 

s Displacement at beginning of braking period 

w Velocity at end of free fall period 

Wo Initial velocity of approach 

C Velocity of exhaust gases of rocket motor 

U Velocity change due to rocket motor 

v Velocity at end of a braking period 

R Radius of moon 

K Gravitational parameter of moon 

h Displacement at end of braking period 

0 Subscript designating any one ON-OFF period 

n Subseript designating last of ON-OFF periods 

m Subscript defined in Appendix II 

Introduction 


A lunar landing is an example of an exclusively 
rocket-powered descent of a space vehicle without 
aerodynamic influence either to assist in the landing 
operation or to hinder in dangerous aerothermodynamic 
heating. If the vehicle is to touch down with zero ve- 

' Received October 1960. 
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locity, then cutoff of the engine must coincide wit 
zero velocity and zero altitude. If the engine is to I 
ignited only once, then the altitude at which ignitio 
must occur is determined by the approach velocity 
This, in turn, requires coincidence of a velocity, a 
altitude and the ignition signal of the engine. Meetin 
these demands of a total of six parameters 1s a sevel 
requirement upon any guidance scheme, particularl 
so if the available thrust is nonvariable. This is obviov 
without an involved mathematical treatment. It 3 
not so obvious, however, that a multistep approae 
relaxes the severity of the scheme to such an exter 
that fairly simple guidance instrumentation become 
feasible. There are several important properties of th 
Multi-Step-Scheme from which to determine its use 
fulness for instrument carrier landings and manne 
landings. Of particular importance are for example 
the ease and economy with which willful lateral ma 
neuvers can be performed, (location guidance), th 
degree of error coupling between lateral and vertica 
guidance, the degree of interference of the rocket flam 
with the guidance sensors and the sensitivity of th 
scheme to the inaccuracies of the guidance sensors 
This report compiles the results of an investigation int 
the progression of errors of a multi-step termina 
guidance scheme based on the following assumptions: 
(1) The thrust is nonvariable. 
(2) The engine is ignited by an altitude sensor 
(3) The engine is cutoff by an inertial velocity 
sensor, Le., after a velocity change has bee 
effected. 
Guidance errors are considered which result from th 


following: 


(1) The altitude, s; , at which the engine is ignited 
is in error by +As, . 

(2) The velocity change, m, which cuts off th 
engine, is in error by +A, . 

(3) The thrust, F,, is uncertain to with +AF; 

(4) The initial velocity of approach, wo, is un 

certain to within +Awp. 

The equation of motion for a braking period canno 
be handled easily because the gravitational componen 
is a nonlinear function of altitude and the mass of th 
vehicle is a function of time. If the gravitation is as 
sumed to be constant the equation becomes linear an 
can be treated explicitly. This approach is taken in th 
analysis. The error introduced thereby is accounte: 
for by consideration of the gravitational acceleratior 
g, of the celestial body (luna) as an “iIndependen 


riable”, with the uncertainty, Ag. This then is a 
fifth error. 

- The effect of these five errors on the final velocity, 
%, , expresses itself in Av, and is quite severe. It places 
very stringent requirements on the permissible mag- 
nitude of the initial five errors. 

_ Earlier work on the subject, (1), (2), indicated that 
an immediate relief can be gained by allowing a small 
free fall period between cutoff and touchdown. 

If a certain impact velocity, w; , is allowed, the effect 
of the aforementioned five errors on w,; is decisively 
less than on v,. The individual errors are reduced by the 
mere insertion of this free fall period; this is shown 
mathematically in Appendix II. Of course, a price 
must be paid. This is the increase of velocity from 1, 
tO UW, . 

One might go one step further and reignite the engine 
after the free fall period in order to eliminate or drasti- 
cally reduce the velocity, w:. This would permit a 
ow touchdown velocity to be linked to the advantage 
f an error reduction afforded by the introduction of 
he free fall period. There is, however, no change in 
he progression of errors. v2 1s Just as strongly affected 
y the initial errors as w,. See Appendix II. Another 
eduction of initial errors could be gained by addition of 
second free fall period, resulting in a final velocity, 
». This would again reduce all errors accumulated 
p to the second free fall period. Therefore, we reason 
hat if this ON-OFF sequence were continued for a 
arge number of times all errors could be drastically 
educed if not practically eliminated. The question to 
» answered, then, is how the errors progress through 
he chain of ON-OFF periods. In other words, if there 
re n ON-OFF periods, it is necessary to determine 
vhat effect an error introduced during the zth preceding 
eriod would have on the nth period. The equations 
f motion used for investigating this error progression 
re derived in Appendix I and are listed as equations 
(12) through (20). 

The present investigation considers a constant 
hrust level during each period but allows for different 
hrust levels from period to period. The lunar gravita- 
jon is approximated during the powered phases only, 
aot during the free fall. The error incurred in approxi- 
ation of the lunar gravitation is accounted for by its 
ntroduction as an independent error in the error 


During each period four errors are introduced. Two 
stem from wrong measurements and two from wrong 
ssumptions. Measurements are taken of the altitude, 
' and the velocity change, uv. The errors are As and 
Au. Assumptions are made in thrust, /, and in lunar 
zravitational acceleration, g. The errors are AF and 
Ag. All periods are alike except for the very first one 
which has an additional, fifth, error in the initial ve- 
locity of approach, wo . This error is Aw» . 


The velocity at the end of the nth period, w, , then 
is a function of: 


(1) m altitudes plus their errors (s; + As; » Sot Ass, 
* 5 8n ae As) 


(2) n velocity changes plus their errors (u + Au, 
Uz + Aug, -+* , Un + Aun) 


(3) » thrust values plus their errors (7, + AF, , 
Fy + AF, , a , Fe AF) 


(4) n gravity values plus their errors (g; -+ Aq, , 
gz + Age, -**, gn + Agn) 


(5) one initial velocity plus error (wy + Aw). 
For the sake of abbreviated writing we shall use: 


(s + As) meaning (s + As, % + As,---, sy 
+ As,) 


(u + Aw) meaning (um + Aw, ue + Au, +++ , Un 
ar Attn) 


(¥ + AF) meaning (Ff, + AF,, F, + AF,, --- 
1b ae AF,,) 


d 


(g + Ag) meaning (g + Agi, go + Age, °°: , 
Grea NGa)e 


Developing w, in a Taylor series with first order 
terms only, as shown in equation (21) of Appendix IT, 
gives: 

Wnls + As, u + Au, F + AF, g + Ag, wy + Aw) 

= w,(s, u, F, g, wo) + (As)A 

+ (Au)C + (AF)D + (Ag) E + (Aw) B. 
Hquations (22) thru (26) of Appendix II show what the 


five A-terms stand for namely: 


(As)A = (As) Ai, + (As2)Amn + -+- + (As,) Ann 
n—1 
= os As; A in + CAS: Are 
t=1 


and analogous for the other coefficients except for equa- 
tion (26) 


( Aw) B = ( Awo) Bin. 
In equation (32) the Taylor series for w, thus becomes: 


Wn(s + As, uw -+ Au, F + AF, g + Ag, wo + Aw) 
n—1 n—1 


= wn(s, u, F, g, mo) + Dd, (Asi)Ain + DD (Au) Cn 
1=1 w=1 


mil n—1 


+ >) (AF) Dn + 2, (AG;) in 4 (Aw) Br, 


pS 


4. CASi Ann + Ate) Cre sF CAF) De ae (Agr) Eve 


The A-terms in equations (32) are easily interpreted 
as to their physical meaning. The last four terms, with 
only the index n, represent the errors introduced during 
the last period itself. The term with coefficient B is 
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representative of the error present in the last period 
because of an error in the initial velocity of approach. 
Finally, the five terms with the >> represent the errors 
that are carried down into the last period from the 
errors introduced during all the preceding periods. 

The coefficients A, C, D, E, B, with their respective 
indices can be resolved to: 


a Pe: ee )e 
Wn Wn-A 


Ge, = Ci a (1 = o ) ’ 
Wn Wn-1 


Din = Dud (1 ~  )a 
Wn Wn-1 

Bin = Pa (1 = ) & 
Wn Wn-1 
1 Un 

By, = Bu — (1 a ) = 
Wn Wn-1 


»= wT (“*), m =1,2,---(n—¢—1). 


(See equations 69 through 74 of Appendix II.) 

All of these coefficients are of the same mathematical 
structure. In other words, all errors progress alike 
through the ON-OFF periods. 

In the equations (69) through (73) a choice has been 
made with respect to the constant value of the gravita- 
tional acceleration of the moon during the braking 
phases. It is advantageous to choose the constant 
g to be the acceleration existing at the moment of cut- 
off of the rocket motor. The reasons are explained in 
Appendix II. The braking phases are thus always cal- 
culated with too large a gravitational acceleration 
meaning that the ambiguous + sign with Ag is re- 
solved into a definite minus. By use of equations (69) 
through (73), the Taylor series of equation (32) be- 
comes as shown in equation (75), which follows: 


w,(s + As,u+ Au, F + AF,g + Ag, wo + Awy) 
= w,(s, u, F,g, Wo) + L ¢ = a 
Wn Wn-1 


[ oAs;A “i co) -|- peers ro) + > AR.Ds Pp 
= Dig. a p aR AwBy & | ++ (As, ) Ana 
ta CM Can aie (AT) Dae = (Ag) Ban 


e=wIl(t — m=1,2,---,(m—i—1) 


Wi+m-l 


© = w,; n—-i~—-l=0. 
Equation (74), above, defines the @-terms. The nine 
error terms of equation (75) are examined critically 
in the following section. 

Equation (75) readily allows the computation of 
the so-called partials or error coefficients of a given 
descent program. These partials are weighing factors 
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which relate the deviation of the final touchdoy 
velocity from the nominal to the deviation of the 
variables from nominal. As an indication of the valuag 
information that they supply; note the following th ' 
partials for the firing altitudes in the descent Progra 
shown in Figure 1. The lower the period the mor 
critical the partials become. There is nearly a two orde 


| 
of magnitude difference between them. 


O19 — 0,050-100* al 


Os) v 

ae -35.10| 

OS. m 

saan NCIS: a 
Os, m : 


Critical Evaluation ; 

In a purely mathematical sense all errors exces 
the ones occurring during the last period could b 
eliminated if u,/Wn+ = 1. (See equation 75). How 
ever, this can be fulfilled only if w,4 is known aceu 
rately, i.e., by measurement of w,; . This renders th 
condition of this equation trivial in a physical sens 
since, if w,1 1s known at a certain altitude, previou 
errors lose significance. 

However, a guidance scheme that measures wy, 
with its associated altitude while neglecting the meas 
urement accuracies in the upper periods might b 
adopted. In that case, the question remains whethe 
the errors accumulated in the earlier braking phase 
would allow close enough approach to the surface s 
that measurements for the last period could be carrie 
out safely and so that the remaining altitude would b 
neither too high nor too low for the actual landing t 
be made. This physical statement interpreted i 
mathematical terms of equation (75) considers wha 
conditions minimize the individual terms in the bracke! 

Consider the series of products &, common to al 
terms, and the individual coefficients A ,; , Cw, DF 
Hii, By. & could be very effective in reducing th 


8: = 278 km -— wo = 2500 m/s 


T uy = 1800 m/s, F = 90,000 Newton 


hy = 158 km -|- ») = 803 m/s 


1 


so = 15.8km -\- w, = 1031 m/s 


t ue = 1050 m/s, F = 90,000 Newton 
he = 2125 m —|- » = 21 m/s 
83 = 226m -|-w. = 814 m/s 


ij us = 84 m/s, F = 40,000 Newton 


hs = 69 m -|- v3 = 1.3 m/s 


0 


w3 = 15 m/s 


Fig. 1. Descent program. This is an arbitrary program whic 
is not optimized or selected with any special conditions j 
mind. It is for illustration purposes only. 


Ititude of errors. It becomes small if the same 
ndition as that mentioned above 


A. 


; met for all or any one of the other periods, that is if 


Uitm — il 
Wi+m-1 


but even if this were only approached rather roughly 
he effect would be very pronounced since the product 
eries of relatively small factors decreases rapidly. 
‘he effectiveness of ® decreases with the lower periods 
ince the product series becomes shorter. Many ON- 
JF periods would, therefore, be advantageous for 
educing errors. For the period immediately preceding 
he last one the effectiveness of © is lost and the co- 
ficients A;;, Ci:, Dix, Hi: must be considered for 
ninimizing the errors introduced during this period. 
his latter point will not be discussed in more detail 
ince all the necessary conclusion can be drawn from 
he subsequent discussion of the coefficients in the 
ast ON-OFF period by simply changing the indices. 

The errors introduced during the last period are 
epresented by the last four terms in equation (75). 
Yne may dismiss the very last one, (Ag,,) Hn», since 
he last period will be so close to the surface that the 
unar gravitational acceleration may safely be assumed 
-constant. This means Ag, = 0. 

The first coefficient, A,, (equation 94), presents a 
lear cut situation. It cannot go to zero. This means 
hat the altitude error when the rocket motor is fired 
or the last time is there to stay. A,,, can only be mini- 
aized by a large w, . This, however, is severely limited 
yy the demanded softness of the landing. Hence, the 
reatest efforts must be made in accurately measuring 
Ititude before firing the engine for the last time. 

Next is C,, (equation 95). It could go to zero if 
» = 0. This is expressed in equation (101) in terms 
f the velocity change, namely: 


ae 0, a a a, a + On tn 
n 
large w, would also be beneficial. 

The last term is D,, (equation 96). Strictly mathe- 

iatically speaking it could go to zero, but for all 
ractical reasons it cannot since the condition for 
), = 0 is only met by a ratio of u,/e = 0 (103, Ap- 
endix II). A compromise w,/c will have to be chosen. 
large thrust to weight ratio would generally be 
worable although at the same time one would prefer 
nw exhaust velocity (notice ¢’) which means low 
pecific impulse. The dependency of D,, on ¢ is rela- 
vely weak, however. 
This evaluation does not yet account for the measure- 
ent of w,_1 which in turn introduces an error (AwWn_1). 
Ihe coefficient with (Aw,1) would be B,, and is 
| 


j 


analogous to By . 


Mp as 
3, ae cas eee EOS 
Wn 
For w,1 = u, the coefficient B,,, is zero. In other 


words the error introduced by measuring the velocity 
Wr-1 before the.last braking starts is not critical. 


Conclusions 


From the standpoint of minimizing the effects of 
imaccuracies of measuring devices the following phi- 
losophy of a lunar terminal guidance with constant 
thrust engines deserves attention. The vertical descent 
should be broken up into a number of ON-OFF periods. 
The guidance effort should be concentrated upon highly 
accurate instrumentation for the last period, mini- 
mizing the errors by proper choice of the velocity 
change and fully exploiting the allowances in the landing 
velocity. For all other periods, the errors can be held 
down with rather simple instrumentation by proper 
choice of the velocity changes, the number of ON-OFF 
periods and the spacing of the ON-phases to the OFF- 
phases. 

The guidance for the last period demands highly 
accurate instrumentation. In particular it requires an 
altimeter for determining the altitude of firing, a 
velocity meter for measuring the velocity before firing 
and a velocity meter for measuring the change of 
velocity once the rocket motor burns. 

This latter velocity meter is certain to be of the 
inertial type and is readily available. The other velocity 
meter and altimeter must be selected with due con- 
sideration for the possible dust stirred up from the 
surface by the rocket flame. Since the last period will 
certainly be initiated below 1000 meters altitude a 
doppler type velocity meter and a radio altimeter 
would be familiar electronic instruments which would 
have none of the problems associated with high alti- 
tude altimeters and doppler equipment. 

If the allowed touchdown velocity is high enough 
(in the order of 15 m/s) the velocity measurement 
before firing can be omitted altogether, thus requiring 
only an altimeter and the integrating gyro. 

For the upper periods a crude altimeter would suffice 
for initiating the braking phases. A horizon sensor 
which is required for establishing the local vertical 
possibly could supply the necessary altitude signals. 
The velocity data for the upper periods may be taken 
from calculated trajectory data and integrating gyros 
thus requiring little instrumentation. 

The thrust level in the upper periods is constant 
during each period but there are no paramount. re- 
quirements for any specific thrust level from the stand- 
point of guidance errors, although a high thrust to 
weight ratio is always advantageous (compare equation 
96). There is also freedom to change the thrust level 
from period to period. The foregoing analysis allowed 
for a different thrust level in each period. This particular 
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point will probably have little meaning for landing 
instrument carriers but a manned vehicle is certain to 
be restricted in the permissible deceleration and a 
reduction in thrust from higher to lower periods for 
the purpose of keeping the thrust to weight ratio within 
narrow limits may be attractive. 

The vertical guidance would impose little restric- 
tions on the lateral maneuvers during the upper periods. 
The lateral maneuvers would enjoy a large degree of 
freedom within the limits of fuel economy. 

The unavoidable errors of the last period must be 
kept within bounds by the length of the final free fall 
before touchdown; in other words, the velocity of 
landing, w,, is the principal means of holding the 
errors down. For any one of the guidance errors the 
velocity dispersion at touchdown is always inversely 
proportional to w, . 

The most important are the altitude error and the 
thrust error. The error in velocity change and the one 
introduced by measuring w,—1 are of second order sig- 
nificance. Theoretically they can be eliminated alto- 
gether whereas the first two cannot. 

If the demands for the softness of the landing exceed 
the minimum landing velocity necessary for keeping 
the errors small (compared to the vertical displacement 
and the vertical velocity of the vehicle), the constant 
thrust engine cannot be maintained. The alternate 
solution would be to aim with the scheme described 
so far (having constant thrust in each braking phase) 
at a touchdown at around 200 meters above the surface 
and switching then to a second smaller rocket motor 
with variable thrust; the thrust would be controlled 
by velocity and altitude continuously. In this case the 
instrumentation is especially concerned with the ioniza- 
tion effects of the flame and the dust problem. Approxi- 
mately a 10:1 range in thrust variability will be re- 
quired. However, it is doubtful whether it is worth- 
while to attempt an extremely soft landing (landing 
velocity in the order of 1 m/s or less in contrast to 10 
to 20 m/s for an ordinary soft landing) without a firm 
knowledge of the lunar surface. Measurements of the 
altimeter and the velocity meter would have to be so 
accurate that only very specific ground properties 
would allow us to obtain these accuracies. Present day 
knowledge of the lunar surface definitely does not aive 
us the required information as to the degree of rough- 
ness, dustlayer coverage, flatness, solidity, ete. 
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Appendix I 


Equations of Motion 


Assumptions: 
(1) The celestial body to which the vehicle descen¢ 
has no atmosphere (luna). 
(2) The engine is first ignited at an altitude whie 
is small compared to the radius of the celestial bod 
(luna). bs 
Because of assumption (2) the surface under th 
vehicle is considered infinitely large and the motion i 
vertical (y) and horizontal direction may be treate 
independently. Only the vertical motion will be deal 
with. 


F 
M = M(t) = M, — M(t) 
g = gly) 
~j=g-—* 
M, — Mt 
+in Mo gm M : ; 


g is a function of altitude y. This makes it hard 
work with equation (3). We will assume that g is con 
stant during the braking phase and account for th 
inaccuracy by considering g as a ‘“‘variable”’ in the erro 
analysis. The assumed error, Ag, will be determined 
by the difference between the actual g and the constan 
g value chosen for the braking phase. 

Vertical velocity, y, and altitude, y, as a function 6 

time, ¢: 


eS a 


y = w + gt — ein! (4 
M, = Mt : 
Y¥ =S$ — Wt — 95 
(5 
$s of 4 ( ~ ue) In | 
a M, — Mt 
M 
i Clm oe 
M, — Mt (6. 


Vertical velocity, y, and altitude, y, as a function 0 
the velocity change, u, effected by the braking: 


2 (7) 
for 0 < Y < wy 


Bs my 8 (1 — o**) 
: Mic Be, 
en ee (8) 


4 Moe’ eel ery ss aa 
C 


wr the free fall following the braking phase the gravi- 
tional acceleration is not approximated by a constant. 


2 


— ;. R 

ee = 9 ty = goo dy (9) 

. sa ox[ ttl 4e lor Wem y <1 10) 
y y if Y 


immarizing the essential equations that express the 
rtical motion in terms of the velocity change, w, and 
troducing the general index 7 one finds for one brak- 
« period plus one subsequent free fall: 


Se, Wi 
ui] 
hi Vv; 
. | 
| 
i 
aw; 
if —u;/c ¢ 
2; = i ad — (1 = " ) =U (12) 
M i€ —tu;/¢ 
h; = 8s; — wi a Gd — ¢“!*) 
= ea re (13) 
Mic —u;/e Uji —u;/e 
ee eset lt a er 1 ey Sy a 
am F. ale eee, = | 
i; = 2K ( : E402 Vi. (14) 
aah Toe 


d the same set of equations for the very last one of 
» periods, n, and the next to last period, n — 1, 


Sn—-1 Wn—2 
Un-1 
Pea 
Ana =e Un—-1 
| 
| 
| 
| 
| 
| 
Sn 7- Wn-1 
‘i 
| 
Un 
F,, 
hy, ran Un 
Sn41 am) 
M nz € 
Un—1 = Wn—2 a: gna a7 = 
n—1 C15) 
(1 = ee 2) — Un-1 
M n—1) © one oO 
hen Oh Ue) ae ve (1 — ‘aur 
2 2 
gn—1 Moin) a Wn AON ) 
og eager Se 1 = @ n—1 
nae ( ) (16) 
Mo =) c =u / Urea 
n 1 e Un— 1/C — nm—1 MAKE 
cs ei \ ) Cc : 
2 sal sil il 2 
Vi, ay = IK — (Nae 
1 E | + 1 in) 
M, nC ame c 
Vn = Wr-1 = Yn (wh a : ) — Un (18) 
Pr 
M n€ Un, /C 
hy, SS perl r, (1 ad nt ) 
_ Gn Mine C ty, /O\ 2 / 
ez Nell reas elas 
| oe oer 
Mone’ ja as nt) = Un el 
ith Cc 
Wn = 2K a =| = Vn . (20) 
Sn+1 Sie 
Appendix II 
Error Progression 
Abbreviations: 
s+ As = &+ Asi, Soot Sahat 5, Susi eS 
au Au = um + Au, Us + Aue , Un a- An 
F+ AF = F,+ AP,, Fy + AF, , , FP, + AF, 


gtAgHntn, Gu, eAGs , dae a On 
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Taylor series for w, with first order terms only: 
Wrls + As, u+ Au, F + AF, g + Ag, wo + Aw) 
= w,(s, u, F,g, wo) + (As)A + (Au)C (21) 
+ (AF)D + (Ag)E + (Awy) B. 


The last five terms in equation (21) are: 


(As)A = (As) Ai, + (Asy) Aon + 
= (22) 
att (Asp) Ann = E (49,) A. = Ase ane 
(Au)C = (Ath) Cin, + (Ata) Con 
n=l (23) 
Ar CSU) ton ad [s (Au) | =F Ain) Onn 
(AF)D as (AF\) Din =F (AF 2) Dan at 
nt (24) 
a (APT) Din = |= (AF )Die| or CATED Das 
(Ag) E = (Agi) Lin + + (Ago) Hon + 
i (25) 
+ (390) Bae = |Z (4g:)Bon | + (Spe) Ba 
(Aw))B a (Awy) Bin. (26) 


The capital letter symbols are the error coefficients, 


also referred to as “‘partials”’. 
Ao a Ave = oe (27) 
at Na rate) 
pe a eee a (30) 
By, = Tr (31) 


€ 


Thus the series of equation (21) becomes: 


Wal s + As, w+ Au, F + AF, g + Ag, wy + Aw) 


= Wns, U, Ei, qY, Wo) 


n—1 n—1 n—1 
ae 2 (As;)A im =F x (Au;) Con = 2 (AF;)D,,, 


n ‘a 99 
+ D | (Agi) Lin Re) 
++ (Aw) Byn + (AB) Airc sr (Aas, IC; nn 
+ (AF, ) Die + (ONG eee 


To obtain the error 
through (20) 


coefficients differentiate (15) 


Vn On IK On» OWn-1 


ei SS ‘ 99 
Wn OS; U7 Ne (33) 


A in — ane a 9 Are ae 
Wn Nr? OWn_-1 OS; 
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(ae Wn seghttl =n 
Wy OUg Walle? OWye OU; 
Dee Vn Phy Wn ixn ( 
si Wn OF, Wnln2 OWn-A OF; 
Em a Un OVn, K Oh, OWn—-1 in ( 
Wn 0g: Wal? OWn-1 0g: 
B, = 22 On Ke One Giles ( 
Wn OW) | Wan? OWn1 IWo 
Oy On, OW OWa at te." (38 
OS; -% OWn1 OS; OS; 
OUn Wn, OWn—-1 es OWn-1 (3¢ 
Ou; OWna OU; du; \ 
Oy, = OVn OWn-1 _ OWA ( ( 
OF, = OWn—-1 OF; OF; . 
On Un, OW 1 ist OWn—1 (4) 
Ode OW = 0G; 0g; 
Wn «Yn OWna  OWn a ( E 
JW)  OWn_1 OW dw : 
KGa Gas ( 
Wale? 2 
: 
(38) through (42) into (33) through (37) and obser 
ing (43) : 
OWA Vy see Oh, | 
A in = 
OSs; ca ut Wn OWn-A 
OWn-1 Un n 
Cn 
Ou; [i +2 Wn OWn— aie 
~ n— n Sed Oh, 
DS ee ak + g | (44 
Wn Wn OW 1 
oh 
Be oe yn — —l1 Vy i n | 
Sanat te ee Yn OWn1 
By = OWy— 1 v n h, 
Jwo a aa 1 
dh, My Cc Uy/ v + Un soe ne | 
= =_>- — aad ag 1 — ginl® a gee n n—1 4 
ae F. ( ev") ; (45 


(45) into (44) 


IWn =] da; [ 9 , 


A in, 
Os 7 Wn ie Gn 


(Wra — 


* 
ibn) =e Un (1 — In ) 
Gn 


; OWaat aes * 
Ce = 7; 3 = (Wnt = Un) Se Vn 1 a LE 


OU; Wh 


OW, zt Ga * 
Din = OF, = in ie (Wy—1 = tia) = Vy (1 a a) (46 


OWrarilig * 
Be eat ie n Wn— i nl ) i Jn 
0g; Wy ie Jn ( n—1 la =- l n 1 


OWn-1 | 
By, a 


* 

Gn ( Gn 
Wn-1 — U Dil Ve — re 

Oy ae oe n) + 


= 


| Equations (46) provide the arguments for choosing 
e value of g, . It is obviously advantageous to select 


In = Ca (47) 
gi = gi*. (48) 


Equations (46) then reduce to the following set of 
uations for the error coefficients: 


ape t =! en (1 —— ) (49) 
- Os i Wn Wn-A1 
B (0 ee (i ~ =) (50) 
3 Ou; Wn Wn-1 

OW, Wn Un ) 

Din = ae cere ol 
is OF; Wn ( Wn-1 \ 
a eee at Ove (1 - =) (52) 
; 0g: Wy Wn-1 

jpn eap ae (1 aa ). (53) 

OWp Wn Wr-1 


Thus dw, has been reduced to a function of Owns. 
1 the same way 0W,_1 can be reduced to a function of 
Ww,» and so on until the 7th period itself is reached. 
nm other words, from: 


A na = OW, — OWT Wn-1 (1 = ae (54) 


OSs; Os; Wn Wn-1 
| 
ve infer 
4 \= OWn—1 a OW,—2 Wn—2 (1 had fay (55) 
£44(n—1 ) 
ee OS; OS; Wn—1 Wn-2 


r in general connotation 


ae GW ay Wee U(n—k f 
i) = OW, k = (n—k—1) “ (n—k—1) i (n—k) , (56) 
OS; OS; W(n—k) W(n—k—-1 


his is also true for Cin_py , Dian—~ and Hin. This 
s true for all integers k except for: 


a = (57) 


»— k = 71 isa special case requiring special attention. 
‘or the time being we stipulate as ‘‘standard coeffi- 
ients’’: 


+ ed (58) 
OS; 
Ow; 
Sige ing (59) 
Cai OU; 
Ow; 
=a (60) 
Di or; 
ee (61) 
0g: 
OW 


Notice that the last equation does not have the index 7. 
For the sake of ease of the mathematical description 
of what is to follow it is expedient to rewrite the equa- 
tions for the coefficients using the index 7 instead of n. 
In other words the rest of the text will be more easily 
explained by looking forward from the ith period to 
the nth period instead of going backwards from the 
last, the nth, period to the 7th period. In this connota- 
tion then the coefficients for the periods following the 
ath period are (taking only one coefficient as an 
example) : 


Aitiuy = Ax ua (1 = faa (63) 
Witt Wi 
W; Uj 
Ae = Aig = 1 — te), 64 
L(i+-2) (i+1) Wiss ( Witt ( ) 


In general connotation this is: 


Witm-1 Uitm ke 
Alen =A Asn) Sanaa (1 Br aa ag od) he 
Witm Witm-1 ( 65 ) 
mM = 1,2. ane 


This is the equivalent to (56). Combining (63) through 
(65): 


Wi Wign \ Witt Wit2 
A itm) = Ai 1 ae Le 1 es ae 
Wi Wi J Wi+2 Witt 


(66) 
ae W(i+m—-1) (i oy U(i+m) ) 
W(i+m) Witm—1) 
Ww; Witt Wit2 
Ashes) = Ai * Loe oe —-— 
Ww it+m Wi Ww +1 
(67) 
Ui é : 
. (1 _ —) m = 1,2,--- (mn —72). 
Witm—1 
And for the nth period finally evolves: 
Wi Wi Ure 
Awe ieee (1 = es 
Wn Wi Witt 
(68) 


(1 -— ) 
Wn-1 


Resultingly the error coefficients in equation (32) 


become: 
Ala Aler i (1 = fs) @p (69) 
Wn Wr-1 
if Un 5 
On = Cy — ¢ —. )e (70) 
Wn Ww n—L 
Date Dae (: - a) & (71) 
Wn Wn-1 
7 1 Un >) 
Ein = Hy, — 1 a ® (72) 
Wn, Wr-1 
1 Un © 
Bae — By =a ( aa ) QD, (73) 
n Wn-1 


N 
wa 
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The operator ® being: 


Wi--m— 
(n —i—1) (74) 


forn —7—1 = 0. 


m=1,2,-°°, 


© = W; 


Enter with (69) through (73) into (32) again and 


OUT ee eer ( 
OWo 
* 
tee = dns ( 
Wn 
C = Wn 0,1 = gn*/Gn) _ Cake Un (¢ : 
nn > Un Wn In Wn 
Doce Wn, On on Dna a Un) + mel = On*/2Gn) 
OT OR Wn 


(oo 
obtain the final form of the series of equation a) 4 Gat Use hy) 
with which the analysis was started: Wn F'n 
w,(s + As, u + Au, F + AF, g + Ag, Wo + Awo) = wrls, u, F, g, wo) 
n—1 n—l . 
hue, (1 a ae b As; Ay ®+ 2S Maa be os AF, Dis + om Ag; Big ® + Aw) Bu *| (7 
Wn Wn— 
== (As, ) Aga ae (Au, ) Cay = (AF,,) Die = (Agn) Bon 
The last four terms in equation (75) contain the error i 2 Wn In*/2Gn(Wna — Un) + Un(1 = gn/2Gn) (95 
coefficients of the very last period only. These are: s On Wy 
K — 0u1 gi*/gi(wo — m) +2101 — n/gi*) , 
= B (93 
Ann aay Wy, h,2 (76) ” ~~ ow Wo Wy 4 
ga in On, a 1S fla (77) Introducing g, = gn* as well as equations (81) throug 
> i OUp eet, Na OU (84) and using with (91) also equation (18) and th 
4 following expression: 
je Dn On i is , Oh», (78) 
Wyn OF, Wy h,2 OF, 2g Myce A 
Danse ay hy) = Un ae Ue 2c 
SO: a Koh, (79) F 
nn Wn Dn Wn hie O9gn -[Un — (Wp—1 = Un ] _ (Wn—r — Un 
Bu “i Vy Ov K ; Oh, ; (80) . 
Wi OW) Wi hy Aw equations (89) through (93) evolve into the fin 


Equations (76) thru (80) can be resolved as follows: 


0 a Gn M n€ 
u ol (Vn + Un — Wnt) = : 
Jy DS 
M oe 
~~ Yn Lon, ane 
—l= a a e 1 
OVn eels ; _ Gnt wD ‘ 
OF, a 7 (On ate Un — Waa) = Fie a (82) 
On a Z (Dn =n Un — Wn-1 ) = bn ( 83 ) 
O0n In 
Ov . 
eS oe | 2 
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expressions for the error coefficients of the very las 
period: 


SAS = Hs Ce = ee | 
& w, , w, (94), (95 
1 Jn Mon C Un —u,/c 
Dae = Wit io 
dee. re Wn | c : # (96 
Bion pea re age : [Wn ry we ie ae | (97 
29n Wn ‘ ‘ 
oy aed 
(om = Se RES sat (98 
Wy 


Kquations (94) through (97) also hold for Ass, 
D;, and H;; by proper replacement of indices. 
Conclusions from (95) and (96): 
= is symbolic for logic KQUIVALENCE 
~ is symbolic for logic NEGATION 
Vv is symbolic for logic DISJUNCTION 


(Ge ra 0) = (Un = 0) v (Wn =  ) (99 


Wn = © is impractical 


Il 


a 
3 

= 

| 
2 
Il 


lie 
(Din = 0) = [e 
(6) 


[un = Wht = Jn tn] 


E = Wr 


eee 


oy | (101) 


= Wr-1 35 In tak 


(102) 


ene 0| = fe = 0) (103) 


=f Thus for all practical reasons: 


~ (Ds = Oe f (104) 
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Technical Notes 


Comments on “Anisotropy Gls 
Escape Velocity from the 
Moon” 


Howard S. London’ 


In a recent article appearing in this journal (1), Gold 
concludes that a difference of 30% exists between the escape 
velocities on the near and far sides of the Moon due to the 
gravitational influence of the Earth. Calculations based on 
Jacobi’s integral of the restricted three-body problem do 
not support this conclusion; an anisotropy does exist but it 
is no more than one-sixth of one percent. 

Mr. Gold compares the minimum velocity required at 
the near surface of the Moon to reach the libration point 
lying on the Earth-Moon axis on the near side of the Moon 
(a libration point being defined as a point where a body 
will remain in equilibrium if initially placed there at rest) 
with the velocity required on the far surface of the Moon to 
escape to zero velocity at infinity. Such a comparison is 
meaningless, since a libration pomt can be reached equally 
well from either side of the Moon. The proper comparison 
therefore is between the minimum velocities required to 
reach this point from the near and far surfaces, respectively. 
These velocities are easily calculated from Jacobi’s well- 
known integral of the restricted three-body problem (2). 
Assuming that the mass of the Earth is 81.45 times the mass 
of the Moon, the diameter of the Moon is 3478 km, and 
that the Earth and Moon are separated by 384,400 km, it is 
found that the required velocities on the near and far surfaces 
of the Moon are identical to six significant figures (= 2.31946 
km //see). 

Even if one chooses instead to compare the velocity re- 
quired on the near surface of the Moon to reach the libration 
point which is on the near side, and the velocity required on 
the far surface to reach the libration point which is on the 
far side, the difference is quite small. These velocities are: 


Viregs surface = 2.31946 km/see 
Viertsuntacce ="erocol2 km/sec 


Thus, even on this basis of comparison the anisotropy is 
only about one-sixth of one percent rather than 30% as 
estimated by Gold. 

Furthermore, the approximate formula Eq. (2b) of (1), 
derived by Gold in (8), gives Vnear surface = 1.889 km/sec, 
which is 18.5% lower than the value of 2.319 km/sec as 
determined from Jacobi’s integral. This discrepancy is due 
to the fact that the analysis in (8) is based upon a model of 
the Earth-Moon system in which the Earth and Moon are 
stationary; such a model is inadequate since it neglects the 
Coriolis and centrifugal forces due to the rotation of the 
Earth-Moon axis. 
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Geodetic Sub-Latitude and | 
Altitude of a Space Vehicle’ 


Robert H. Gersten”’ ° 


This technical note outlines a procedure for determin 
the geodetic sub-latitude and altitude of a space vehicle fro 
its geocentric position. The geocentric and geodetic su 
latitudes of the vehicle, denoted by ¢’ and @ respective 
are related by the conventional equations (see Reference 
1 and 2 at end of article): 


(C + H) cos = r cos g’ 


(S + H) sin g = rsin @’, ( 
where, utilizing a., e, and f to represent the Earth’s equ 
torial radius, eccentricity, and flattening and r to represe 
the separation between the vehicle and geocenter, 


a es oF i 


(1 — e sin? ¢)!? 


e= ap | 


and where H is the height above the geoid expressed in 
same units as a, (see Fig. 1). Multiplying Equation (1) a 


— Space Vehicl 


North Pole 


~— =& —> 
Geocenter 


Fig. 1. Meridonal cross section of the Earth showing tt 
vehicle’s geocentric distance 7 and geocentric and geodet 
sub-latitudes $’ and ¢. 


Equation (2) by sin ¢ and cos ¢ respectively, and subtractir 
yields: 


(C — 8) sin g cos g = r sin (¢ — 9’). (: 


' Received October 1960. 
* Senior Scientist, Astrodynamics Department. 
* Norair Division, Northrop Corporation, Hawthorne, Cali 


1 Equation (3) 


a, e 
(@ — e sin? )12 : 


C-S=eC= 


ming a binomial expansion on the above 
4 C — § = ae*(1 + fe? sin’g) 


sre powers of e greater than the fourth are truncated. 
stitution of the above into Equation (4) yields: 


sin (6 — ¢') = e sin 9 cos A(1 + de? sin’g) (5) 


Le 


y 


A, e 


° 


a 

R, 
order to maintain all terms through the e* terms in Equa- 
1 (5), it is only necessary to find expressions for sin ¢ and 
@ as a function of ¢’ and e, which are expanded through 
é term (since the right side of the Equation is multiplied 


ae Alp : ; ee 
a . This is accomplished in the manner indicated 
! ? 


OW. 
Substituting Equation (5) into 


cos (¢ — ¢’) = [1 — sin? (@ — g’)}'? 
| expanding, yields 
cos (¢ — g’) = 1 — 3 sin’ + $e? sin‘g. (6) 


stitution of Equations (5 and 6) into the trigonometric 
ntities: 


@ = sin [g’ + (¢ — 9’)] 

= sin g’ cos (¢@ — ¢’) + cos @’ sin (¢ — Q’) 
cos [¢’ + (@ — 9’)] 
cos g’ cos (¢ — g’) — sin g’ sin (g — ¢’) 


s 
I 


ll 


sin ¢ sin ¢’ + cos @’(e sin ¢ cos @) 


cos 6 = cos @’ — sin g’(e sin G Cos @). 


previously explained, it is unnecessary to keep terms 
taining e raised to higher than the second power (or e€ 
ed to higher than the first power) in the above in order to 


maintain accuracy through the e! term. Substitution of the 
expressions for sin @ and cos ¢ above into the sin ¢ and cos 
terms on the right side of the above yields: 


sin g = sin g/(1 + e cos?¢’) (7a) 

cos g = cos ¢/(1 — e sin?9’). (7b) 
Now substitution of Equations (7) into Equation (5) yields: 
sin (¢ — 9g’) = esin 9’ cos g’[1 + « — (2e — 4e*) sin’g’] 
where 


Ae @ 
(ee 
rr 


The above may be expressed in terms of flattening by: 


sin (¢ — ¢’) = aE sin 29’ + f? sin 49’ (* — i) | ; 
r r 
Similarly, eliminating C and S from Equations (1, 2, and 
3), the geodetic altitude, H, may be obtained as 
H =r — a,({l1 — $e?(1 + e) sin?’ + Je2(e — +e?) sin‘¢’] 
where again 


Ae 
€ = 


r 


or, in terms of the Earth’s flattening, 


: I? 2. (pn PAN 
H =r—a|1-—f sin? ¢’ —=sin? 29’ [ — —— , 
| ik ey (° ‘) | 


See pages 166 and 167 of (2). 
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Letters to the Editor 


Comments on ‘“The Influence of Heart Action and 
the Circulation of Blood on Manned Satellite 
Attitude Control,” 


Gentlemen: October 28, 1960 


In the Fall 1960 (Vol. VII, No. 3) issue of the Journal of 
the Astronautical Sciences the paper by Lawrence R. Zeitlin 
entitled “The Influence of Heart Action and the Circulation 
of the Blood on Manned Satellite Attitude Control” contains 
a, fundamental error in mechanics that leads to a false con- 
clusion about the requirements of the satellite attitude con- 
trol system. 

Specifically, Mr. Zeitlin first calculates the effective angular 
momentum of the blood circuit of a man. He then calculates 
the angular momentum of the satellite for any angular 
velocity. He then, erroneously, equates these two angular 
momenta and solves for the angular rate of the satellite. 
Furthermore, he then states, erroneously, “Over a long period 
of time, a significant amount of energy would be required 
from the attitude control system to counteract this torque”, 
produced by the “reaction to the angular momentum of the 
astronaut’s circulatory system”’. 

Assuming that the satellite and the enclosed astronaut 
have been initially stabilized, the only time the satellite will 
change its angular momentum (start rotating) in an inertial 
reference frame is when the angular momentum of the 
astronaut’s blood changes in the inertial reference frame, and 
then the change in the satellite momentum will be such as to 
keep the angular momentum of the satellite plus the angular 
momentum of the astronaut’s blood constant in the inertial 
reference frame. If the astronaut is kept relatively motion- 
less in the satellite, as postulated by Mr. Zeitlin, the only 
justification for equating the two angular momenta is to 
assume that the astronaut’s blood has stopped circulating. 
Under this condition the satellite will indeed begin to rotate 
(with angular momentum equal to the angular momentum 
of the astronaut’s blood before it stopped circulating), how- 
ever, an angular momentum wmpulse from the attitude 
control system can stop this rotation. No further energy 
expenditure from the control system will be required to keep 
the satellite from starting to rotate again. 

Considering a more general situation (also more palatable 
to the astronaut), assume that the astronaut is free to move 
about inside the satellite. He is thus able to change the diree- 
tion (in the inertial reference frame) of his blood’s angular 
momentum. The satellite will rotate in the inertial reference 
frame so as to maintain constant total angular momentum 
of the system (satellite plus astronaut). Note that angular 
momentum is a vector quantity and that a constant vector 
has constant magnitude and direction. If the astronaut moves 
from one orientation to another the satellite may start to 
rotate but an angular momentum impulse from the attitude 
control system can stop this rotation, and it will not start 
again until the astronaut moves again. If the astronaut moves 
around considerably the total energy involved in all the 
momenta impulses may become significant if a non-con- 
servative attitude control mechanism is used, such as re- 
action gas jets. A highly suitable conservative attitude 
control mechanism is available however, in the form of re- 
action flywheels. These reaction wheels provide a third com- 
ponent of angular momentum in the overall system that can 
balance the astronaut’s blood momentum changes and thus 
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prevent the satellite from rotating when the astronaut change; 
orientation in the satellite. 
Henry M. Musat, JR. 
Senior Engineer, Research Laborator 
Bendix Systems Division, 
493 East 31st. Street, 
Chicago 16, Illinois ~ 


November 30, 19€ 
Gentlemen: 


Mr. Musal is perfectly correct in stating that an error was 
made in the original paper. This was an error in computatic 
however, rather than a “fundamental error in mechanies. 
We simply did not use the same assumptions regarding init 
stabilization of the satellite. The error was corrected in 4 
revised copy of the paper, which unfortunately, did 
arrive until after the publication deadline. The correcteé 
material should read: 


of blood per minute; or about 7200 Kg. in a 24 hour day 
(1). At any time, he has the equivalent of a fluid flywhee 
within his body with a mass of approximately 5 Kg., s 
radius of 5 cm., and rotating once a minute. This flywhe 
has a small, but finite amount of angular momentum 
which, in accordance with the principle of conservation 6 
momentum, would exert a rotating torque on his bodi 
in the opposite direction. Since the blood flows clockwise 
looking from the front, the torque tending to rotate th 
body would be counterclockwise. : 

What would the effect of this rotating torque be on a 
small manned satellite? ; 

The angular momentum of a body is defined as th 
product of its rotary inertia and the angular velocity 
around a given axis. For a flywheel with its weight con 
centrated in the rim, the moment of inertia (J) = MR 
where J is the mass and R is the radius. Hence the angula 
momentum of the equivalent blood circuit is: 


(5 X 10%) (25) (6.28) 
60 


= MRW = 


= 1.31 X 10' gm em? radians/second. 
For a 1000 Kg. satellite, of approximately the shape 
of the Mercury vehicle, the moment of inertial along its 


axis is that of a right cone of altitude H radius of base R, 
and mass J/7. That is: 


3 
I=M—R 
10 


For the Mercury vehicle, 


M = 1000 Kg 
R = 100°cm 
H = 300 em 


Hence: 
I = 108 X 104(.3) = 3 & 109 gm em?. 


According to the law of conservation of momentum, 
the angular momentums of the two rotating masses must 


dl 


4 equal. Hence: 


1.31 X 104 gm cm? rad./sec. 3 X 10° gm em? W 


Ang. momentum of blood Ang. momentum of vehicle 


_ W = 4.36 X 10-6 rad./sec. 
: or 
2.5 X 104 deg./sec. 


In a 24 hour day, a 1000 Kg. satellite would rotate 21° 
40’ about its axis simply in reaction to the angular mo- 
mentum of the astronaut’s circulatory system. 

I fully agree that, given a state of initial stabilization, the 
nly time the satellite will change its angular momentum in 
n inertial reference frame is when the angular momentum 
f the astronaut’s blood also changes in the inertial reference 
ame. Mr. Musal gives a backhanded justification for equat- 
i momentii of the satellite and astronaut by stating that 
the astronaut’s blood stopped circulating the satellite if 
abilized would begin to rotate with an angular momentum 
jual to that of the angular momentum of the astronaut’s 
lood. Being of a more sanguine nature, I preferred to start 
ith a satellite-astronaut system of zero initial angular 
omentum and determine the maximum rotation of the 


system attributable to the angular momentum of the as- 
tronaut’s blood. No assumption of initial stabilization was 
made. If uncorrected, the rotation proposed by Mr. Musal 
as a function of the death of the astronaut would be equal 
to the rotation computed above, but in the opposite direction. 
It’s just another view of the elephant. 

As for the more general case, the angular momentum of 
the astronaut’s blood changes from moment to moment. 
States of sleep, excitement, digestion, deep thinking, and 
activity all produce changes in blood pressure and volume 
of circulation which can be reflected in changes of angular 
momentum. Try taking your pulse after resting, eating, or 
chasing your secretary around the desk. Pulse (and angular 
momentum) may vary as much as 3:1 under these conditions. 
The attitude control system, whether conservative or non- 
conservative, would have to correct continually for these 
angular momentum changes if precise attitude control is 
required since these changes are not of a discrete nature. 

By the way, if the astronaut can orient his seat in any 
direction, he can serve as his own reaction wheel. 

LawRENcE R. Zerrin 
Dunlap and Associates, Inc. 
355 W. Olive Ave. 
Sunnyvale, California 


Book Review 


Review of “Interplanetary Communications”, Chapter 2 
F Advances in Space Science, Vol. 1, Academic Press, New 
ork and London, 1959. 

This is an excellent survey of the subject by Dr. J. R. 
ierce and Mr. C. C. Cutler of Bell Telephone Laboratories. 
; is natural that the Bell Telephone System would explore 
ve technical and economic feasibility of using space to 
ve its long distance communication problems such as the 
‘ansatlantic telephone problem. The highlights of a high- 
rade study of low and high altitude communication satellites 
gether with pertinent charts and formulas are most wel- 
yme. 

Besides the fine treatment of the use of space vehicles to 
crease the number of channels and reliability while greatly 
ducing the cost of Earth communications, this Chapter 
so considers the problems of interplanetary and interstellar 


communication. The treatment is from the modern sophisti- 
cated information theoretic and practical point of view. 
Modulation methods are compared with the ideal, but as 
yet unrealized, technical optimum. Technical-economic 
optima for choice of transmission frequency are given. It is 
interesting to see that the granularity difficulties which would 
arise for communication at X-ray frequencies are not serious 
for a wavelength of 1 cm. Here the authors estimate 80 
quanta per cycle bandwidth. 

This summary of ideas on space communication, with its 
sprinkling of profound remarks on what is important, is 
highly recommended to both the expert and the layman. 


Joun L. BarNnus 
Systems Corporation of America 
1007 Broxton, Los Angeles 24, California 
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Format of Technical Papers 
4 for AAS Journal 


The Editors will appreciate the cooperation of 
hors in using the following directions for the prepara- 
n of Manuscripts. These directions have been com- 
led with a view toward eliminating unnecessary cor- 
spondence, avoiding the return of papers for changes, 
4 reducing the charges made for “author’s correc- 
fanuscripts 

Papers should be submitted in original typewriting 
f possible) on one side only of white paper sheets, and 
iould be double or triple spaced with wide margins. 
Owever, good quality reproduced copies (e.g. multi- 
th) are acceptable. An additional copy of the paper 
ill facilitate review. 

ompany Reports 

The paper should not be merely a company report. 
‘such a report is to be used as the basis for the paper, 
propriate changes should be made in the title page. 
ists of figures, tables of contents, and distribution 
sts should all be deleted. 

‘itles 

The title should be brief, but express adequately the 
ibject of the paper. A footnote reference to the title 
10uld indicate any meeting at which the paper has 
een presented. The name and initials of the author 
10uld be written as he prefers; all titles and degrees or 
onors will be omitted. The name of the organization 
ith which the author is associated should be given in a 
sparate line to follow his name. 


.bstracts 

An abstract should be provided, preceding the intro- 
uction, covering contents of the paper. It should not 
sceed 200 words. 


leadings 
The paper can be divided into principal sections as 
opropriate. Headings or paragraphs are not numbered. 


llustrations 

Drawings should be made with black India ink on 
hite paper or tracing cloth, and should be at least 
puble the desired size of the cut. Each figure number 
ould be marked with soft pencil in the margin or on 
ie back of the drawing. The width of the lines of such 
rawings and the size of the lettering must allow for the 
scessary reduction. Reproducible glossy photographs 
ce acceptable. However, drawings which are unsuitable 
r reproduction will be returned to the author for re- 
rawing. Legends accompanying the drawings should 
2 typewritten on a separate sheet, properly identified. 


ecurity Clearance 

Authors are responsible for the security clearance by 
1 appropriate agency of the material contained in the 
ypers. 

fathematical Work 

As far as possible, formulas should be typewritten. 
reek letters and other symbols not available on the 
pewriter should be carefully inserted in ink. Each 
ich symbol should be identified unambiguously the 
st time it appears. The distinction between capital 
1d lower-case letters should be clearly shown. Avoid 
mfusion between zero (0) and the letter O; between 
e numeral (1), the letter 1, and the prime (’); between 
pha and a, kappa and k, mu and wu, nu and », eta and n. 
The level of subscripts, exponents, subscripts to sub- 
ripts, and exponents in exponents should be clearly 
dicated. 


Greek Alphabet 

A a_ alpha (a) N vp nu (n) 
B £6 beta (b) = & xi (x) 
T y gamma (g) O o  omicron (6) 
A 6 delta (d) Il x pi (p) 
E ¢« epsilon (&) P p_ rho (r) 
Z ¢ zeta (z) > os sigma (s) 
H 7» eta (8) T > - tan (t) 
6 6 theta (th) ae upsilon (u) 
1% ota (i) & ¢¢ phi (ph) 
K «x kappa (k) X x. chi (ch) 
A 2 lambda (1) VW vy opsi (ps) 
My mu (m) Q w omega (5) 


e (epsilon) = strain 
g,s (sigma) = stress 
7 (tau) = shear stress 
p (mu) = micro 


p» (mu) = Poisson’s Ratio 


Complicated exponents and subscripts should be 
avoided when possible to represent by a special symbol. 

Fractions in the body of the text and fractions occur- 
ring in the numerators or denominators of fractions 
should be written with the solidus. Thus: 


cos (7x/2b) 
cos (7ra/2b) 


is the preferred usage. 

The intended grouping of handwritten formulas can 
be made clear by slight variations in spacing, but this 
procedure is not acceptable in printed formulas. To 
avoid misunderstanding, the order of symbols should 
therefore be carefully considered. Thus: 


(a + bx) cos t is preferable to cos t (a + bx) 


In handwritten formulas the size of parentheses, 
brackets and braces can vary more widely than in 
print. Particular attention should therefore be paid to 
the proper use of braces, brackets, and parentheses 
(which should be used in this order). Thus: 


{[a + (b + cx)"] cos ky}? 
is required rather than ((a + (b + cx)") cos ky)?. 


Equations are numbered and referred to in text as 
(15). 
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